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1.0 SUMMARY 


Engineering preliminary design methods for approximating and predicting the effects of 
viscous or equivalent viscous-type damping treatments on the free and forced vibration of 
lightly damped aircraft-type structures are developed. Similar developments are presented 
for dynamic hysteresis-viscoelastio-type damping treatments. It is shown by both engineering 
analysis and numerical illustrations that the intermodal coupling of the undamped modes 
arising from the introduction of damping may be neglected in applying these preliminary 
design methods, except when dissimilar modes of these lightly damped, complex aircraft-type 
structures have identical or nearly identical natural frequencies. In such cases it is shown 
that a relatively simple, additional interaction calculation between pairs of modes exhibiting 
this “modal resonance” phenomenon suffices in the prediction of interacting modal damping 
fractions. The accuracy of the methods is shown to be very good to excellent, depending on 
the normal natural frequency separation of the system modes, thereby permitting a relatively 
simple preliminary design approach. This approach is shown to be a natural precursor 
to elaborate finite element, digital computer design computations in evaluating the type, 
quantity and location of damping treatments. It is expected that in many instances these 
simplified computations will supplant the more elaborate ones. 



2.0 NOMENCLATURE 


2.1 Notation 

A cross-sectional area of rod, in 2 ; an arbitrary constant 

a coefficient in interaction quartic equation; membrane dimension, inch 

B coefficient in resolvent cubic equation; an arbitrary constant 

b coefficient in interaction quartic equation; membrane dimension, inch 

C coefficient in resolvent cubic equation; damping coefficient, lb-sec/in 

c coefficient in interaction quartic equation; damping coefficient, 

lb-sec/in 

D plate flexural rigidity, lb-in 2 ; coefficient in resolvent cubic equation 

d coefficient in interaction quartic equation; damping coefficient, 

lb-sec/in 

E Young’s modulus of elasticity, lb/in 2 

F function symbol; force 

/ force, lb 

G shear modulus of elasticity, lb/in 2 

I area moment of inertia, in 4 ; mass moment of inertia, lb-sec 2 -in 

i ordinal number; a subscript 

J torsional section constant, in 4 

j ordinal number; complex operator, v^T 

K stiffness, spring rate, lb/in 

k stiffness, spring rate, lb/in 
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mass, lb-sec 2 /in 

mass, lb-sec 2 /in; an ordinal number; a subscript 

an ordinal number; aspect ratio 

an ordinal number; a subscript 

naught, a subscript 

generalized coordinate 

a response quantity ; a generalized force 

a generalized coordinate 

a response quantity; a parameter 

an ordinal number; a subscript; a ratio 

an ordinal number; a subscript 

membrane tension, lb/in 

time, sec 

a displacement, inch 
a displacement, inch 
a displacement, inch 
a Cartesian coordinate, inch 
a Cartesian coordinate, inch 
a Cartesian coordinate, inch , 

an ordinal number; a subscript; a constant 
an ordinal number; a subscript; a constant 
a parameter 
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A dilatation 

6 a damping coefficient, lb-sec/in; perturbation symbol 

C fraction of critical damping 

7] loss factor 

6 torsional displacement, radian 

A, A characteristic number 

fj, mass per unit length, lb-sec 2 /in 2 ; mass per unit area, lb-sec 2 /in 3 

v frequency ratio 

£ axial coordinate position 

p mass density, lb-sec 2 /in 4 ; frequency ratio 

cr damping per unit length, lb-sec/in 2 

Q forcing frequency, rad/sec 

u> frequency, rad/sec 

2.2 Symbols 

dot, differentiation with respect to time 
bar, amplitude of 

* asterisk, complex conjugate 

V 2 del squared, the Laplacian operator 

arrow, vector quantity 
tilde, a modified quantity 
f ( ) integral, the integral of ( ) 

| | the magnitude of, the determinant of 
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square matrix 


[ ] 

U row vector 

{ } column vector 

d ( ) differential of ( ) 

d ( ) partial differential of ( ) 

div ( ) the divergence of ( ) 

grad ( ) the gradient of ( ) 



3.0 INTRODUCTION 


Aircraft, spacecraft, and especially rotorcraft airframes are subject to steady forced 
vibrations due to a variety of rotating or oscillatory type mechanical and aerodynamic 
systems. These steady forced vibrations can become severe when resonant or near resonant 
conditions occur in the airframe. For example, in the case of rotorcraft, excitation frequencies 
include the rotational frequencies of the main rotor and tail rotor (both different) and 
harmonic excitations at integer multiples of their blade number times the fundamental 
frequency. In the case of two-bladed rotors, for example, excitation frequencies at rates of 
once per revolution, twice per revolution, four times per revolution, etc. are commonplace. 
In the case of propeller /rotor-type systems, similar families of excitations exist, but the 
difficulties can be compounded if the propeller/rotor operates at one rate of revolution in 
hovering flight and another one in forward flight. 

Of crucial importance in the case of resonant or near resonant forced vibrations is the 
fraction of critical damping associated with the particular mode which is responding. At 
resonance there is a direct, inverse relationship between the magnitude of the response and 
the magnitude of the fraction of critical damping of the mode. If damping can be increased 
by an order of magnitude, then the response is reduced by an order of magnitude, etc. In 
the case of airframes, the inherent damping levels are small and of the order of a few percent 
of critical damping or less. Accordingly, augmenting the normally small levels of inherent 
damping can be very beneficial when a severe steady vibratory response is attributable to a 
resonant or near resonant condition. 

Since weight, cost, complexity, etc. are among the primary concerns in airframes, then 
whenever a “fix” or corrective addition of damping is indicated, a simple but accurate 
preliminary design type of engineering method of analysis is to be desired. This method 
should be capable of both rapid and relatively simple, but nevertheless accurate engineering 
predictions of modal damping. This is especially important in guiding the structural 
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designer-dynamicist in determining the location, type, and quantity of energy dissipation- 
damping treatments when the airframe is sustaining severe forced vibrations. 

The objective of this report is to develop a simple, preliminary design-type of analysis 
and methodology which can accurately predict the modal damping associated with damping 
treatment. It is expected that very detailed, lengthy, and complex finite element type 
computations may also be performed, especially in the development of a new airframe 
structural design. However, the preliminary design method developed and presented in 
this report is intended to be a precursor to such finite element type computations. The 
designer can quickly determine the efficacy of localized damping treatments to within 
acceptable engineering accuracy prior to undertaking a much more detailed and complex 
design computation. It is also expected that the accuracy of this preliminary design method 
will frequently obviate the need for any other computations. 
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4.0 ANALYSIS 


4.1 The Single Degree of Freedom Oscillator With Damping 

The special case of a dynamic system whose vibrations can be adequately represented by 
a single degree of freedom is considered first. This will permit the definition and development 
of various damping concepts in a simple manner before dealing with the general case of a 
damped structural dynamic system with many degrees of freedom. 

Consider first the case of viscous damping where the damping force is proportional to 
velocity, but opposite in sense. The governing scalar differential equation of motion is 
mx+cx+kx = fei^, where the complex exponential represents a simple harmonic excitation 
of amplitude / and frequency fi. The parameters m, c, and k represent the system mass, 
damping and spring constants, respectively. 

In the transient case when / = 0, the general solution takes the form x(t) = xe^. A 
substitution of this leads to the characteristic polynomial A 2 4- (A) A + (jj^j = 0. Defining 
the natural frequency of the oscillation as w, u 2 = , and its fraction of critical damping 

by C = ( 2% /ra fc ) ’ characteristic equation becomes A 2 + 2£uA + ui 2 = 0. This equation 
yields the characteristic number A and its complex conjugate A* as follows: 

A = -v + 3\ A -C 2 ) . A* = -ui ( C - j \[l - C 2 ) , 3 2 = -1* 

The transient solution then takes the form 

x(t) = e-^ (a COS U) \fl-c 2 t + B sin a; V^t), 

where A and B are constants to be determined from the initial state of the transient 
oscillation, where 

x(o) — x 0 and x(o) = v 0 . 

In the case of steady forced vibration the transient oscillation decays rapidly and after several 
cycles the steady state response is given by x a .s.(t) = Xs.s.e^ ■ Substitution above results 
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in the equation 


x s.s. 


f 


(k — rafi 2 ) + j (cQ) 
Introducing the previous definitions /, fi, and £ gives 


*s.s. 


(i 

) e-M 

1 

ED 

j] 

+ 

R (g)f 


where 


(j> ~ tan 


-1 


2 CI 

( 8 ) 


1-1 

(n> 

r. 


Defining the static displacement x s . s . static = (//&)> the steady state forced response is given 

by 


x s .s.{t ) = 


X S.S.static e ‘ 




> - (S) 2 ] + R(g)f 

Now consider the case of dynamic hysteresis, where the damping force is proportional to 
displacement, but opposite in sense to the velocity. Employing the imaginary operator j as 
a “phase shifter” or “differentiator,” the governing equation can be written as follows: 


mx + k( 1 + jr})x = feP^, 


where 77 , the “loss factor” will be seen later to be a measure of energy dissipation by dynamic 
hysteresis. In the transient case when 7 = 0, the general solution can be written as x — xe^, 
so that the characteristic equation takes the form 

A 2 +u/ 2 (l +jri) = 0. 

Accordingly 

A = ±ju(l+ jri) 1 / 2 . 

Employing DeMoivre’s theorem, 


( 1 + ^) 1/2 = (1 + 7 ? 2 } 1 / VW 2 ) 
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where 


6 = tan 1 77, 

A = -w(l + t? 2 ) 1 / 4 ^sin ? ± j cos 0 . 

In the typical engineering case of 77 < 1 and, generally, 77 -C 1, 

A -(*)’)• 

Thus the loss factor divided by two may be interpreted as equivalent viscous damping with 
^Equivalent = r //2- As an alternative to the foregoing transient vibration analysis employing 
the operator j, the governing equation with dynamic hysteresis may also be written as 


.. fkrj\ . , 

mx + I — I x + kx = 0. 

\u J 

Taking the solution once again in the form x(t) = xe^, the characteristic equation follows 
as 

A 2 -f- ujt)\ ■+■ uP — 0. 

It is seen at once that CEquivaient = r ?/2- 

In the case of steady-state forced vibration with dynamic hysteresis employing either the 
formulation 

mx + k( 1 + jrf)x = /e^ 4 , 
or 

mx + f— ^ x + kx~ /e^ 4 

leads to a steady-state solution with comparable results. First with the operator j: 


x s . s .(t) = Xg - S ' s ^ tlc - 




v / [ i -(§)] 2 +» 2 


where 


1 

(f> — tan 


V 


m 
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The equivalent viscous damping then follows as 


^Equivalent 


’ll 2 


m/u)\ ’ 

so that strictly speaking, only at resonance (Q/u) = 1 is CEquivalent — ??/2 (however this is 
the case of greatest practical interest); in the case of the other formulation 

, I V (§) 


x s .s.(t) 


X S.S. Static e ' 


h(8) 2 ] + MS)] : 


(j) = tan 


-1 


Ms)' 


The equivalent viscous damping then follows as CEquivalent = r?/2. 

That r] is a measure of energy dissipated per cycle can be seen as follows: the dynamic 
hysteresis damping force is given by 


/d 


ampmg 


= —krjx 


( 6 )- 


Computing the work done against /damping P er cycle of simple harmonic response at 
amplitude x Q and frequency fi, 

(2tt / ft) (27r/f2) (2 tt/0) 

w/cycle = — J fddx = — J fd^dt = kr)x^l J cos 2 = wkrjXg. 

0 0 0 

By contrast when the damping force is modeled as viscous with /damping = — Equivalent^ 
the energy dissipated per cycle follows in a similar manner as 

w/ cycle — 7rcQr 2 . 

Equating the two expressions for work per cycle 

__ kri 

Equivalent ' q ' 


Dividing through by m 


Equivalent 


m 


2CEquivalent^ - ( m ) (q) ~ ( Q ) * 
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Then 


CEquivalent 



This is seen to be the same as the result obtained above when forced vibration with the 
operator j is taken as the basis of defining the relationship between the loss factor 77 and 
CEquivalent- ^ f° be noted that it is primarily at or near resonance that CEquivalent fakes on 
special significance. Accordingly, from an engineering perspective the relationship is taken 
as a universal one that 


CEquivalent 


V 

2 


Various other forms of energy dissipation are also of interest during forced vibration. 
Once again special interest centers on the response at or near resonance. Also, the cases 
of relatively small amounts of damping are assumed, so that despite the nonlinear aspects 
and character of the dissipative force, the steady state response exists, is dynamically stable 
and for practical purposes is simple harmonic in time. Employing the concept once again 
of energy dissipated per cycle as a basis for deducing and defining the equivalent viscous 
damping coefficient, the following summary can be made for the most frequently encountered 
cases of damping: 

(1) Dynamic Hysteresis 


^Equivalent — > CEquivalent 


( 2 ) Coulomb Friction 


_ (?) (*£“) , 

Equivalent — q > s Equivalent 


L(S) 


(V 

\ ir > 


( 8 ) 


(3) Velocity Squared Damping (/Damping = (jf|)) 

^Equivalent = (3^) (P*o)> CEquivalent = [(3^) (“) (™)] 5 
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(4) Velocity Cubed Damping (/damping = -ax 3 ) 

3 ^ 2 2 > r/3 \ (tfxlw 

Equivalent — x oi ^Equivalent ~ I \8 / \Vmfc/ J * 


4.2 Free Vibration of Systems With Viscous or Equivalent Viscous Damping 

The governing matrix differential equation of free motion of systems with viscous or 
equivalent viscous damping is 

M;{*} + [ c]{x } + [k]{x} = {0}, 

where {x(t)} is the vector displacement, [m], [k], and [c] are, respectively, the system mass, 
stiffness, and damping matrices. The [m] and [k] matrices are symmetric and positive 
definite; the [c] matrix is symmetric and either positive definite or positive semi-definite. 
Consider a solution in the form 


«*)} = 

where and A j are, respectively, the complex modal vector and the associated complex 

characteristic number, where i = 1,2, ... ,1V, and 

{x}i = (x} iReal + i{®}i ImaginMy , j 2 = -i, 

’ ^*Real $ ^Imaginary * 

There are also N complex conjugate vectors {x}| and their associated characteristic numbers 
A* which satisfy the differential equation. The assumed form of solution leads to a system 
of algebraic equations which follows in matrix form as 

[a |H + Aj[c] + [fc]j {x}i = {0}, i - 1, 2, . . . , N. 
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Pre-multiplying by the transposed complex conjugate modal vector |xj* yields N scalar 
equations 

A? (|®|* M (®}i) + \ (ifl* W {*}i) + (ifli W {*}i) =0, * = 1, 2, • • • , N. 

This is re-written as 

Af + +ujf = 0, 

where 

2 _ ( ^ / 1 \ ( ^ 

These N scalar characteristic equations yield the N characteristic numbers and their N 
complex conjugates which are as follows: 

A* = ~u>i i j ^ » i = 1? 2, . . . , N. 

Orthogonahty relationships which are useful in analyzing forced vibrations are now 
developed. Writing the r th and s th matrix algebraic relations: 

[A?[m] + A r [c] + [&]] {x} r = {0}, 

[AjH + A s [c] + [k ] j {x} s = {0}. 

Converting both of these to scalar relationships by pre-multiplying by [xj s and |x| r , 
respectively, yields: 

A? (|x| s [m]{z} r ) + A r (|^s[e]{x} r ) + (ifjsWMr) = 0, 

A? (l£lrH{x} s ) + As ([x)r[e]{x} s ) + (|x| r [fc]{x}s) = 0. 

Substracting the s th equation from the r th and noting the symmetry of [m], [c], and [A;] 
yields 

( - Aj) ([x] r [m]{z} s ) + (A r - As) ([£|r[e]{x}s) = 0. 
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Factoring out (A r — A s ) when r ^ s, the orthogonality relationship follows as 

(|£jr[c]{£} s ) + (A r + As) (|x| r [m]{x} s j = 0. 

It is to be noted in passing that if A s = A*, then the previously derived result for ( r 
follows from this orthogonality relationship. 

An alternative orthogonality relationship can be deduced by adding the r th and s 1 ^ scalar 
equations above. Noting the symmetry of the [m], [e], and [k] matrices once again, there 
results the scalar equation 

(a? + A^) (|4H{*}s) + (Ar + A s ) (|£)r[c]{£}s) + 2 (|x]r [&]{£}«) = 0. 

Noting from above that 

(Ar + As) (|£)r[c]{x}s) = -(Ar + As) 2 Qx]r[m]{x}s) 
and substituting yields 

[(a 2 + A^) - (Ar + As) 2 (|i|r[m]{i} s ) + 2 (|£jr[fc]{>}s) = 0. 

Simplification results in the alternate orthogonality relationship 

(|fjr [&]{*}«) - A r A s ([x| r [m]{x}s) =0, r ± s. 

It is to be noted in passing that if A s = A*, then the previously derived result for follows 
from this form of the orthogonahty relation. 

In summary when r ^ s, 

([®|r [c]{*}«) + (A r + As) (|i|r[m]{i} s ) = 0 

and 

(|*jr [*]{*}«) ~ A r A s (|£|rM{£}s) = 0. 
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When the s th mode corresponds to the complex conjugate of the r th mode, then 

2£ r qj r = 

\lf]rH{®}r 

and 

(j 2_ f 

\l£fN{®}r > 

It is to be noted that, except in the case of differing modes with closely matched natural 
frequencies, the undamped modes may be employed in calculating u r and, in turn, £r- This 
is discussed at length below and in Reference 1. 

4.3 Free Vibration With Two Different Modes Having the Same 
Natural frequency 

To illustrate consider the case when modes “r” and “s” have the same natural frequency, 
although the modal patterns differ. A simple example of this is provided by a rectangular 
membrane of length a and width b with interchanged nodal lines such as the modes w mn (x, y) 
and w nm (x,y), where 

W r (x,y) = W mn (x,y) = sin zsin y, 

W s {x, y) = Wnm(x, y) = sin (^j a; sin V, 


and 


o /7rar\ 2 /ri7r\ 2 j /T\ 

**■- Fr) + (t) Ju) 


UJ, 


nm 


f mr \ 2 /m7r\2 

+ (x) 


rwi \- 

a ) 


©■ 


where T is the membrane tension force per unit length and y is the membrane mass per 
unit area. Equating the foregoing frequency expressions, a frequency match occurs whenever 
a = b, which corresponds to a square membrane. A more general case occurs when u/ TOn == u rs 
and the membrane aspect ratio N = ( a/b ) satisfies the constraint N = 


„ 2 Jl 
m—r 


s 2 — n 2 


V2 , and 


m, n, r, and s are integers. For example when m = 1, n = 2, r = 3, and s = 4 a frequency 
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match occurs for the aspect ratio N = (a/6) = .8165. Other detailed examples will be 
provided later in the section presenting numerical results. 

Now consider an N degree of freedom system in free damped vibration as above where the 
r th modal vector {x} r differs from the modal vector {x} s , but the r th natural frequency 
u r equals the s th natural frequency ui s . Now expand the modal vectors into an undamped 
mode series. That is with undamped modal vector participation factors or series coefficients 
p n , where 

N 

= ^ 1 Pn } ^Undamped ‘ 
n=l 

The system algebraic equations when {x(i)} = {x}e^ then take the form 

N 

[A 2 M + A[c] + [fc]] {x} = ^2p n [A 2 [m] + A[c] + [fc]] {x}n Undamped = {0}. 

n=l 

Noting that in the case of the undamped modes 


|-u> 2 [ra] + [k] ] {x} n = {0}, n = 1, 2, . . . , N, 

N 

53 Pn [(** + w n) H + A W] {S}ra Undamped = {0}. 

n=l 

Also noting the orthogonality of the undamped modes given by 


M H = lf| Undamped H^JjUndamped = °> ^ * 


the matrix algebraic equation simplifies to 



w 2 0 


( Pi ' 

f°l 

A 2 [J] + X[d] + 



* 

- ? 


. 0 


IpnJ 

loj 


where [/] is the identity or unit matrix and 


. ( li!ll ^Undamped ^Undamped ) 

UtAA — 


*IJ 


M 


n 


r 


i-,3 — 1, 2, . . . ,N. 
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Note also the interactive damping fractions Qj where 


Cij 

Also let v = and r { = (g) . Then 


dij. 


v 2 [I] + v 


d\id\2 

^21^22 

dmdN2 


d iN 

d 2N 

d NN 


+ 


0 


0 


'N 


Pi 
P 2 

PN 


> = 



The characteristic determinant can, of course, now be expanded to obtain the characteristic 
roots of the system. However suppose that r a = rp, the case of different undamped modes 
with the same natural frequency. In this case of the degenerate modes “a” and the 
damped modal vectors are linear combinations of “a” and “/?” • That is 


Pa = 7Pc* + 6pp, 


-/ 


Pp = + PP/3, 


or 


(Pa \ = h ^ifpal 
1 p'pj L e p\ \ppy 


Since only “a” and “(3” interact in this special case, then the two algebraic equations for 
the modes “a” and “/?” decouple from the rest of the system due to this degeneracy. This 
results in a relatively simple and informative quartic characteristic equation for these modes. 
That is 

r 

2 


V 


1 0 
0 1 


+ v 


daa d a p 

dpa dpp J 


. 2 
+ r 2 


1 = 

0 iJJV^j W’ 


and 


v* + {daa + dpp)v Z + (2r^ + d aa dpp - d a pdp a ^ u 2 + r 2 {d aa + dpp)v + r% — 0. 


When 


( d aa dpp d a pdp a ) 


daa d a p 

dpa dpp 


= 0 , 
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v A + {daa + dpp)i/ 3 + 2 rji/ 2 -f r 2 (d aa + dpp)v + r% = 0. 

This quartic is now seen to be factorable into the product of the quadratic factors 

(v 2 + rj ) [i/ 2 + (< d aa + dpp)v + rjj = 0. 

The implication of this is that when matching natural frequencies of different modes occur, 

one of these modes is undamped, while the other has the damping of both. More about this 

will be presented in the section on numerical results where it will be seen that many cases of 

engineering interest satisfy the foregoing conditions, or are numerically comparable to this 

when either ^a/3 j s exactly zero or negligibly small compared to 2 r 2 above. 

“/3a “/3/3 

4.4 Employing the Undamped Modes to Determine the Damped Modes 

It has been shown above that the exact undamped natural frequencies and modal damping 
fractions for the various modes of the dynamic system are given by 

2 _ / 1 N / |*|fl<9{*}i\ 

Clearly the efficacy of an engineering approximation depends on accurately approximating 
the damped modal vector |xjj and its complex conjugate |x||. The undamped modal 
vector |x|j w provides the necessary approximation, provided that the system is lightly 
damped and that two dissimilar modes do not have the same or nearly the same natural 
frequencies, the “modal resonance” case. Express the damped modal vector as the 
z tJl undamped modal vector and a perturbation effect due to damping, {x) Ui and 
respectively, where 

{x}i = {x}ui + Cii{dx}i, 

_ / i \ / |£|m[c]{g}m \ 2 ( \ 

\2vuJ % 

In approximating the mode fraction of critical damping fox use, for example, in forced 
vibration calculations, it is clear that if {6x}i is the same order of magnitude as {x} Ui and 

■ " ; -'.W 
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Qi is very small compared to unity, then {x}^ and {x}* will approach {x} Ui and Cn — Q. 
Rayleigh’s theorem for damped linear systems states that the eigenvalue A^ is stationary 
with respect to perturbations in the eigenvector: that is, 

dXi/d{x}i == 0 . 

Employing the perturbation form for {x}i yields 

dAj ^A}/rf{r}uj q 

d{x} Ui + Qid{6x}i [l + C iid{8x}i/ d{x}ui J 

For very small values of {Cm} and absence of a “modal resonance” which can result in a very 
large { 8x}i , the eigenvalue is seen to also be stationary with respect to the undamped 
modal vector {x} Ui , so that 

d\i/d{x}i — 0 . 

In this case 

— ~ u, Ui ii "b 3 \J^- — C^ i 3 = 1* 

The perturbation vector {6x}i is now examined by expanding the modal vector {x} into 
an undamped modal vector series, 

N 

{x(t)} = } 

i- 1 

where the pi(t) are generalized coordinates or undamped modal participation factors. In 
view of the orthogonality of the undamped modes with respect to [m] and [fc], 

N 

Pi T ^Cii^UiPi *b ^UiPi ~ ^ ~ $ij)Pji 

3 = 1 

where is the Kronecker delta with magnitude zero or unity as i ^ j or i = j, respectively. 
Also, C ij is given by 

< ij = (1/2 w Ui ) ([x|ujc]{x} uj / • 


20 



Cleaxly the perturbation in p{ (and in turn the i th damped modal vector) is of the order of 
Cm, unless a “modal resonance” occurs. 

4.5 Forced Vibration of Systems With Viscous Damping 

The system is assumed to be in steady forced vibration under a periodic excitation. This 
periodic excitation has a typical simple harmonic component {a;(t)} = where {/} is 

the vector amplitude and Q is the excitation frequency; the complex exponential represents 
the simple harmonic variation in time. As in the case of free vibration (x(t)} is the system 
displacement vector, [fc] and [m] are symmetric, positive definite matrices representing system 
stiffness and system mass, respectively, [c] is a symmetric matrix which is either positive 
definite or semi-definite; it represents either actual viscous damping or an equivalent viscous 
damping representation for other forms of energy dissipation. Accordingly the governing 
matrix differential equation is 

M{*> + [«]{*} + WM = {/(*)} = {fh ,nt 


The system is now represented in a canonical matrix form by augmenting the system 
displacement vector with the system velocity vector. The augmenting identity equation 
[m]{x} — [m]{x} = {0} yields the following partitioned matrix format for the system forced 
vibration: 


H i M f{*>! M ! [0] fto 


.M j [o] 


[0] I —[ml 




A still more condensed form results by employing {z 
then becomes 




the governing equation 


[A\{z) + {B\{z} = {F}^\ 
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where 


f mi 


[c] ] [m] 


[fc] j 

[0] 

j r- 

o 

1 W-' 

> , [A] = 

_M | to] . 

, [■ B ] = 

JO] I 

-M. 


In the case of homogeneous free vibration there are 2 N damped modes and their associated 
damped characteristic values {z} r and A r , respectively (i.e., there are N complex modes and 
their associated complex characteristic numbers and their N complex conjugates). These 
satisfy the damped modal orthogonality relationships 

[zJ r [A]{z} s = [£j r [i?]{z} s = 0, r ^ s. 

A modal expansion solution is now developed for the steady state forced vibration of the 
system. In the steady state {z(t)} is represented as 

2N 

{*(*)} = q r (t){z} r , 

r = 1 

where the q r (t) are generalized coordinates or “damped modal participation factors.” 
Substituting above and employing the orthogonality relationships, 2 N uncoupled, scalar 
response equations result. These are 

Hrrq(t) + v rr q r {t) — r = 1 , 2 ,... , 2 N, 

where 

Hrr = (l£|r[A]{z} r ) = (jx| r [cj{x}/) + 2A r (|x] r [m]{x} r ) 

v rr = (|fj r[5]{z} r ) = ([fjr [k] (x} r ) - A? ([zjrMIxfy) 

Qr = (ifj r{F}^ S (|x]r{/}) 

In the steady-state the modal participation factors q r (t) have the solution q r (t) = ij r e J . 
Substitution above yields the complex amplitude q r given by the equation 




The complex closed form solution for steady-state forced vibration then follows as 


2 N 

{x(t) } steady— state = ^ ^ v Qr{ u rr + j^f^rr) {%}r- 

r= 1 

Substituting from above and employing the previously defined modal scalar quantities, 
the steady-state response of the damped system is 

st»dy **«. r=1 V ) [(|f| r [c]{i} r ) + 2Ar ((ijrHWr)] / 

where 


Also 


and 


X r = —&r ^Cr + j \J*L — Cr ^ • 

^ 2= / l£^[fe]{^}r 
y[xl*[m]{x} r 

/ 1 \ / \ 
\ 2 u r ) ^[£j*[m]{x} r J 


\ r 

where the asterisk (*) denotes the complex conjugate modal vector. 
The amplitude of the r th generalized coordinate is 


Qr — Qriyrr "b j^(J>rr) 


This can be rewritten as 


Qr ^"Static 


1 + jm 



where the r th mode static displacement is defined as 



In the case of free damped vibration, the r th mode has the Rayleigh-type quotient for the 
r th characteristic exponent 
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Accordingly the ratio ( =-& — ) follows as 

V ” r Static / 


^ Static 

Similarly for the (r + 1) mode 




( fr 


\ Qr -Hlstatic 


;) = H (a 


-1 


Since the damped modes will occur in complex conjugate pairs for all cases of practical 
engineering interest, let the (r + 1) damped mode be the complex conjugate of the 
damped mode and denote this by an asterisk superscript^ Then 


( ) - 

r!_ 7 Y n> il 

\ ^static / 

r ^UJJ 


-1 


Consider the portion of the steady-state response due to the damped mode and its 
conjugate: 


Qr{ x }r + Qr{ x }r ~ ^static i x }r 1 ~ 3 


-1 


+ ^Static 




or 


Qr 


^Static 

1 ~3 

(5)] {*h +5; SttU . 

i- 

■*(£)]{*}? 


1 ~3 


] 


Continuing to simplify 


Qr 


^Static 

i-j 

’ (§)] 

r “1“ 9rgtatic 

1 ~ J (ft). 

{[ 

H 

Q2 V 

(i 

+ i)} 



Noting that A r A* = and (A r + A*) = ~2( r u r , 


^Static 

l- 

3 (ft)] 

■r+Qr, 

"l 

Static 

1-3 

(ft)] {*>* 

{ 

Hs) 1 2 

+ 3 

2 C,(S)] 

1 


1 Subsidence or critically damped modes will not occur in the aerospace type structures of interest. 
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Then the complete closed-form solution for the steady-state response of the system to a 
simple harmonic excitation follows as 

M*)} S teady— state = £ j |l - 

+j {§;)] } { ?rStefe 1 ~ j (§)] {$}r + ^ta t io 1 ~j (3;)] ’ 


where 


_ Qr 

^Static — „ » 

Qr = [f|r{/}, 
2 

z/j'i* = K.<tt — j 

K rr = |x| r [fc]{x} r , 
M rr = |x| r [m]{x} r , 
X r = —u> r jsj\ — Cr ^ > 


* = QJ 

^Static 7/* * 
rr 


( Iff [k]{*}r 

ui r = 1 


1/2 


q rs 

Or = [ff{/}> 

i/* r = K* r - A * 2 M* r , 

K* rr = \$[k\{xy r 

M*r = |ff M{*}*, 

A* = ~Wr ^Cr — — Cr^ > 

|£l* [c]{x} r 


^ 2 o; r ) 


[ff [m]{x} r y ’ ’* \2u} r J \^|xf [m]{x} r 

The foregoing steady-state form of the solution is now rewritten as follows: 


N 


{*W} S te^y- 8 ta te = E eJ(!!, ' W 


r—1 


where 


^Static [" 

~ 3 (§)] {*}* + ?r static f 1 “ 3 (£)] &Yr 


4 -«)T*K(«r 

= tan -1 

2 & (S) ’ 

b-(£) 2 J 

• 


It is to be especially noted that each damped mode (and its conjugate) respond with dynamic 
magnification and phase shift characteristics similar to single degree of freedom damped 
oscillators. The magnification depends on the ratio of the forcing frequency to the natural 
frequency (Q/u; r ) and the modal fraction of critical damping ( r - It is seen that a resonant or 
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near resonant modal response magnitude is crucially dependent on £ r , the fraction of critical 
damping. 


4.5.1 Approximating the Response of Lightly Damped Systems 


Consider the numerators of the response solution. Defining this as A r : 
1*1 r{f} 


Ar -■ 


K rr ApA 


+ 


K* r - A * 2 M, 


‘"'( 8 ) 1 ** 

[-( 8)1 


{x}*, r = 1,2, .. .,21V. 


Neglecting the imaginary parts of {x} r and {x}* for small damping yields {x} r = {x} 

{*}rundamped- Then 




) 


zM 

-(*)*] 


+ 


1 ~3 


»(*) 


m 


Noting that 




(S)-z 




w r ^Cr +j V 1 ~ 


^-" 0 ( Cr+j v^) 


and that 


(av) ~ + = + 2 iC r\J 1 ~ Cr^ + K 


and 


(£) =(-i-2* v / r^)+2« 

Ifk.tf} 


At — 7r 


J^lru. W(^}ru. y 
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where 


lr = 



[1 - (Sr) VW? 

+3 

Cr (£) 

11 

1 - i (Wl - Cr 2 

J 


+ 


[[i+(S)v^+<? 

+ 3 


m 

[ (Cr\/l ""Cr 

J 


For practical purposes, if Cr < 1, the complex parameter 7 r in the r th damped mode response 
is unity. Accordingly the lightly damped system response can be computed, mode-by-mode, 
with a quasi-single degree of freedom response function R r . That is 


N 


7=1 


where 


Rr= f [£|ru,{7} 

\ l^jru. [&] {^}ry. 


-m 


+ [ 2Cr (£)] 


l 3 ' In?. {/} is the scalar energy or work function of the load distribution {/} acting in the r th 
(undamped) mode and |g|ru.[fe]{^}ru. is the r th generalized or modal stiffness scalar. The 
ratio of these two scalars is analogous to the static deflection of a single degree of freedom 

1 fl r U.{/} : ru- „th 


system. That is, 


{/} is a static or very slowly varying one so that 


is a measure of the r tb mode response when the load 


{/(*)> = 0}^ m = {?}, 


or in effect 


( 


n\ 

— <« i. 

u r J 


4.6 Free Vibration of Linear Systems With Dynamic Hysteresis-Viscoelastic 
Damping 

The governing system differential equation in matrix form is 

H{£} + [k]{x} + j[Sk]{x} = {0}, 
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where {x) is global displacement vector for the iV-degree of freedom system, [m] and [k] 
are symmetric, positive-definite matrices. [£fc] is a symmetric matrix which is positive-semi- 
definite. The complex operator j, where j 2 — —1, represents and effectuates the dynamic 
hysteresis damping force character where these forces have the sense of being opposite 
to the system velocity components, but are at the same time proportional to the system 
displacement components. 

The general solution for free vibration displacements is taken in the form (a:(t)} = {x}e xt , 
where {x} is a complex modal vector to be determined, and A is a complex scalar 
characteristic value associated with the modal vector. This reduces the matrix differential 
equation to the following matrix algebraic equation for the mode: 

A IN + [k] +j[6k]\ {x}i = {0}. 

Pre-multiplying this equation by the transposed complex conjugate modal vector |jg||, a 
scalar equation for A* follows below. 

x i (yiM{*}i) + ([£|i[ fc ]{ 5 }i) +3 ([£liW{£}i) = 0. 

Defining the real scalar quantities 



the scalar equation for A j is rewritten as follows: 

Xi + u?(l + jra) = 0. 

Solving for Aj, 

x i = j 2 vj0- + jm), 

and 

k = + jr)i) 1/2 ■ 
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Employing Euler’s theorem 




1 +j»ft = (l + Vi) e^ tan 


»7i) 


and 


Accordingly 


(1 + i»K) 1/2 = (l + Vi) V4 e^ taa * »?)] . 

A i = (l + » 7 ?) ^ |+ sin Q tan -1 77 ^ - j cos tan -1 77 ^ 


Its complex conjugate A| is 


A i - ~ u i ( x + Vi) / |+sin Q tan 1 Vi') + j cos Q tan 1 77 ^ j 
The mode decay factor is the real part of the characteristic value given by 

(l + Vi) / sin Q tan -1 Vi) • 


A *Real ~ W * 


By analogy with viscous damping where the decay factor is A i Real = — Ciyiscous^’ ^ f°^ ows 

Viscous 


that 


( 1 + Vi ) / sin Qtan 1 v^j , 


Real 

Equivalent 

where Vi is the modal analogy of the material loss factor 77 for a simple structural element 
with complex moduli 

E = E{1 + jv) or G = G( 1 + jrj)- 
Generally 77 is much less than unity, so that for 77 •< 1 


Ci viscous Vi/ 2 ‘ 

Equivalent 

In the extreme cases of synthetic rubber materials (specially compounded and impregnated 
with carbon-black particles to enhance damping), the loss factor can approach unity in 
magnitude. Then for rji = 1 


™ Equivalent 
Viscous 


= (2) 1 / 4 sin ^tan 1 1^ = (2) 1 / 4 sin ^) = 0.4550. 
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It is seen that in this extreme case, the linearized loss factor prediction of the equivalent 
viscous damping fraction of critical damping of 0.5000 overestimates the exact value by 
about nine percent (9%). A graph of the modal equivalent viscous damping fraction over 
the range of modal loss factors from zero to unity is shown in Figure L It is seen that 
negligible error occurs over the practical range of interest from % = 0 to rft = 0.5; where at 
rft = 0.5, CiEquivjjent = -243 and the error of the linearized approximation is only about three 

Viscous 

percent (3%). 

Returning to the definition of 7)i, 

\ [fl *[&]{*}; / 

and expanding the complex modal vector into the representation 

{x}i = {£} iundamped + 

where the real, undamped modal vector {^}i Undamped is perturbed by the complex vector 
as Tji increases from zero. Substituting above, the equation for becomes 

([xjiu + 7?f 

Vi ■= “^7 \ 7 “ . 

Except in the case of a modal resonance, the perturbation vector {&r}j is small, generally 
very small and of the order of magnitude of itself or less as in the case of lightly damped 
viscous systems (discussed above). Accordingly the approximation for % employing the 
undamped modes is seen to be a valid engineering approximation for rfo. That is 
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4.7 Pervasive Dynamic Hysteresis-Viscoelastic Damping 

In the special case when the matrix [8 k] is given by [6k] = r)[k] an exact decoupling 
of the damped modes is possible employing the undamped modes. In this case the matrix 
differential equation for free vibration is 


\m]{x) + (1 +jn)[k]{x} - {0}. 

Taking the damped mode solution in the form {x(t)} = {*}z Un damped e ^’ ^ follows that 

[ A iH + (1 + 3V)[k}] {^}iundamped = W' 


Employing the undamped mode orthogonality condition 


Iflnj. H Msu. = |®|ru. l k ) { a '}«u. = °> r # s > 


then A| + u;|(l + jrj) = 0, where is given exactly by Rayleigh’s quotient 

2 = ( \ 

\lllm.H {»}%./ 

It follows as in the general case that 


= “"d 1 + r; 2 ) l ' ,4 sin Qtan 1 7 1 ] = 

*W»y = w i( 1 + ’» 2 ) 1/4c “Q ttm_1 'l) -“i 


for 77 >C 1. In this special case where the loss factor 77 is more a universal structural property 
of the system, rather than a modal one, the equivalent viscous damping ratio of every mode 
is the same, namely 


^ Equivalent 
Viscous 


v 

2" 
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4.8 Forced Vibration of Systems With Dynamic Hysteresis 

Consider a simple harmonic system excitation at frequency 0 (radians/second) repre- 
sented by the complex exponential vector {/(£)} = {/}©#**. The system matrix differential 
equation is 

[m]{x} + [/c]{x} -I- j[6k]{x} = Q}e? nt . 

The steady-state forced vibratory response can then be represented as {x(t)}sg = 
{x}s.s.e^- Substitution above yields the matrix algebraic equation 

[k] - 0 2 [m] +j[6k]j {x} s s = {/}. 


Now consider the damped modal expansion 

N 

r=l 

and rewrite the foregoing equations 


Y Pr [N ~ ^ 2 N +ip]] {x} r = {/}. 

r=l 

Continuing, this equation is equivalent to 


N 

Y,P r [fcl + SW + A rH ~ ( fi2 + X r) N] {®}r = {/}• 

r=l 

Since 

[[k] + j[Sk] + A 2 [m]j {x} r = {0}, 

then it follows that 

N 

- Y,P r ( fi2 + X r) HWr = {/}• 

r= 1 

Pre-multiplying by the transposed complex modal vector |x| s and employing the dynamic 
hysteresis-viscoelastic orthogonality condition 


[xj r [m]{a;} s = 0 r^s 
32 
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" 1 •• : 





(ft 2 + xty PrM' rr = y r{/}, M' rr = ([xjrH{^}r) , 


and 


pr “" (sTO?) (St) ’ Qr = ^ {7} ' 


The r th and s th damped modes can be interrelated as follows: 


[®|s [[fc] + j[6k] + A 2 [m]] {x} r = 0 


and 


[xjr [k] + j[Sk] + A 2 [m]J {»} s = 0. 

Since [fe], [6fc], and [raj are symmetric matrices, subtraction of the second equation from the 
first yields 

(A 2 - A 2 ) ([£j r [m]{x} s ) = 0. 

Accordingly, since A r ^ A s , the orthogonality condition follows as 


Defining 


and noting from above that 


xl r [ra]{a;} s = 0, r ± s. 


_ ( Qr 

Static - M ,^J1 


A 2 = -w 2 (l +jrjr), 


/ Pr \ _ r 

4 

■ 1 

VPrstatic / L-^ 2 + w r(! + JVr) . 

— S 



and the magnitude of this response ratio, “the damped magnification factor with modal 
hysteresis” follows as the magnitude 

1 


Pr 

P r Static 


-(£) : 


+ Vr 
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4.9 Derivative Operator Formulation For Systems With Dynamic 
Hysteresis 

Consider the free vibration of a structural dynamic system with dynamic hysteresis type 
energy dissipation. We seek the equivalent viscous damping fraction of critical damping for 
the i th damped mode of vibration. The governing matrix differential equation for the free 
vibration is taken as follows: 


The incremental viscoelastic type forces have the phase of viscous damping forces, but do 
not increase in magnitude with increase of natural frequency of vibration (hence the factor 

M- 1 ). 

The general solution to the governing equation for the mode of free vibration is taken 
in the form 

{ret*)} = 


Then 



where 

of = Ku = ((sjlWIWi) , 3 

fiKjj — l^lt {^}i ■— {‘ E }iR ea l d" ^Imaginary’ ^{^Ihmaginary 


Defining 


Then 


8K< 


n 


CilJVniivalftnt * 


M U Ui ‘ ******' 


1 / 1 


^Equivalent 2 l MijU? 


(6Ki i). 


Substituting for u;? 


^Equivalent 


1 (SK U \ 

2 \K U )‘ 
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Defining the i th the modal loss factor as 



It is to be noted that for practical purposes of modal loss factors (and internal loss factors 
also) of the order of 0.50 or less, there is negligible difference in this result and the one 
employing a complex stiffness approach for modeling the viscoelastic-dynamic hysteresis 
effects in the damping of structural vibrations. 


4.10 The Global Equations With Proportional Damping 

Consider the case of so-called proportional damping, when the damping matrix [C] is 
expressible as a linear combination of the system global mass and stiffness matrices. That is 

[C] = a[m] + (3[k\. 

Accordingly the governing matrix differential equation for free vibration of the structural 
dynamic system takes the form 

[m]{a;} + [a[m] + /?[&]]{£} + [&]{#} = {0}. 

In the undamped case the characteristic vector {x} u i and the associated undamped 
natural frequency u> u { satisfy the matrix algebraic equation 

[[*] ~ 4n] {x}ui = {0}, 

and the orthogonality relations 

jfjm M {x} u j = lf]m[fc] {x} uj = 0, i^j 
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In the damped case, consider a solution for the damped characteristic solution in the 
form 

{•e}* — 

If this is vahd, then 

[A? H -I- Aj[a[m] + P[k\] + [fc]j {x} ui = {0}. 

Pre-multiplying by the transposed i th undamped characteristic vector results in this 
matrix algebraic equation being reduced to the scalar algebraic equation 

X i (|f|mM{z}m) + >H +^[^]]{^}m) 

+ (ifJuiWMm) = 0 

It can also be inferred that {x}j = {x} u i, since other potential terms in {ir}j involving the 
modes {x}j ^ {x}% will vanish due to the orthogonality of the {x} u i with respect to [m] and 
[k\. Accordingly the foregoing scalar equation is an exact relationship for determining the 
damped mode characteristic value when proportional type damping is present. It follows 
then that in this case 


\ — 0, * — 1) 2, . . . , N, 


where 


a/? = 
w i — 


_ ( Mui W {x}ui 


|*|uiM {*}«», 

CiS (^)( 0+/5 ^) 

For example when proportional viscoelastic or proportional dynamic hysteretic type damping 
is assumed where 

U>i 


then 


A _ 

’h^uivalent 2^°* 
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5.0 APPROXIMATING THE DAMPING FRACTION 
FOR CONTINUOUS LINEAR SYSTEMS 


5.1 The Damped Rod in Axial Vibration 

Figure 2 illustrates a uniform rod which is built-in at one end and free at the other. It is 
embedded in viscoelastic material capable of dissipating energy during vibration. Neglecting 
the stiffness of the damping material compared to that of the rod itself, a one-dimensional 
wave type equation follows which describes the damped, free axial vibrations of the rod 

AEu" = fjLU + aii, 


with the boundary conditions 


«M) = =0. 


Noting that the natural frequencies of the undamped free vibrations are given by 

(t) (?) 


AE 


w n 


, n — 1,2,..., 


and that 


\fi J { m/t ) \m/ 


the governing partial differential equation is separated with the solution 


u{x,t) = u(x)e^ t . 


This leads to the characteristic equation 


from which 


where 


n^n^n + — 0 


— -w n ^Cn i j \J 1 — i 


Cn 


2mu) n ’ 
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and where 

c=ai and m = yJL. 

Collecting the various definitions and substituting above yields 


Cn — 


1 

c 

(2 n — l)7r_ 

^AEfx_ 


n — 1,2,... . 


Figure 3 illustrates the case of the rod once again, but with a damper of rate c placed at 
the previously free end. The governing equation is now the wave equation 


AEu" — fiu, 


with the boundary conditions 


w(o, t) = 0 

AEu'(£,t ) = ~cu(£,t) 

In view of the foregoing boundary conditions and the solution in the form 

u(x,i ) = u(x)e^ t , 

there results the characteristic equation 

tanh 


(M)< 


'JL 

AE 


+ ^ = 0. 


Since the damping fraction as a function of the system parameters is desired, the character- 
istic exponent is transformed as follows to reveal the real part explicitly. Let 


y ^ (M) =^R + i A /> f - -!• 


Defining 


7 


(v^Ae) 


> 0 , 


^ ix AE / 

employing the trigonometric-hyperbolic function identities 


sinh jz = j sin z 
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cosh jz = cos z, 


the transcendental, complex characteristic equation simplifies to two separate equations, 
each of which must vanish. These are 


cos A/ = 0, 


and 


tanh Ar = —7. 


Employing the foregoing definitions 



and since the admissible values of A / are 


( 2 n-l\ 
= (~) 


n — 1)2 ,-,. , 


the natural frequency of the damped oscillation is 

'2n-l\ /7T\ /AE 


w n 


{2n-l\ / tt \ [AE 

~\ 2 )\e)]l )i ’ 


n = 1,2,... . 


This is seen to be identical to the results for the undamped case. Continuing, 




and the transcendental equation for the decay factor and damping fractions for the various 
modes follows as 


. AE 

Inverting the hyperbolic tangent function 
^Reai = ~ UJ n(n 


tanh 6 


y/JIXE 


= “ 0 )Vt tanh ~ 1 ( 


v^A E) * 
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Since lightly damped systems are the ones of primary interest, an engineering approximation 
follows at once by employing the power series expansion of the hyperbolic tangent function 
for small argument z: 


z 3 z 5 


tanh z = z + — + — + - - (z < 1); 
6 5 


taking 


y/JlAE' 


A «Real ** 



\AE\ ( C \ = _ fc \ = _ /C \ 
)) \ld)~ W’ 


and 


ix J WixAe, 


Cn^n = ( — ) , and 
Vm/ 




)tt_ (vM-e) ‘ 


[(2n-l)7rJ \^/fxAE , 

It is noteworthy that a damper of rate c at the otherwise free end of the rod yields twice 
the damping fraction of the pervasive, uniformly distributed damping c = crl previously 
determined. 

The engineering solution for the damping fraction for the case being considered is now 
shown to be identical to that obtainable by neglecting the intermodal coupling due to 
damping, which is seen to be valid for lightly damped systems. The n th undamped mode of 
free axial vibration of the rod is 


u n (x) = sin ( '—y ( j) x, 


n = 1,2, ... . 


The governing, uncoupled differential equation is by analogy with a single degree of freedom 
system 

^ n Effective^ n ^Effective ~ 0, 


and 


\ ^Effective / 
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where 


M 


^Effective 


Accordingly 


and 


J q n(x)ul(x)dx = fij^ sin 2 (|) a: dx 

^Effective = J Q cr(x)ul(x)dx = (^) sin 2 7r = a 

2C n w n = — » Cn^n — ( 1 j 

rn \TO/ 


fj£ TO 


2 

( c "1 

(2n — 1 )tt _ 

WhAe) 


Cn — 

This, of course, is identical to the formal solution of the result above employing an exact 
closed form solution for a lightly damped system. 

To illustrate this point further by a direct numerical comparison, consider the rod with 


a free-end damper designed to yield a fundamental mode damping fraction £1 


Nominal 


= o.io. 


This is a relatively large fraction of critical damping for an engineering structure. Employing 
the exact, closed form transcendental equation solution, 


CW = 0.1014. 

This implies that the exact solution is 1.4 percent greater than that calculated from the 
engineering approximation. This small error decreases to zero as the nominal fraction of 
critical damping desired is decreased. For example for a nominal damping fraction 

= 0-02 and = 0.02006, 

the approximation error is approximately 0.3 percent. 


5.2 Damped Structural Members Other Than Rods 

Since the governing equations for simple, St. Venant type torsion are also one dimensional 
wave type equations, it follows that by analogy the damping fraction results are identical 
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to those of stretching, axial type rod vibrations. For engineering cases which are lightly 
damped and the damping treatment ranges from x = x\ to x = x% the damping fraction 
follows in general as 

K ^ ra EffecHve j 

Effective ) 

where 



/ P t ( x Wni x ) dx 

Jo 

rl _ rx 2 _ 

^Effective = / <r( x ) d U x ) dx = a / 0^(x) dx. 

do d Xl 


^Effective 

r l 

/o 


Similar formulas for biharmonic or flexural systems such as beams, plates, and shells can 
be written by inspection. They are expected to be even more accurate than those for wave 
equation type systems. This follows from the greater separation of natural frequencies and, in 
turn, still weaker intermodal coupling due to damping. For example, a comparison of natural 
frequency spacing for a simply supported, uniform beam shows a mode number squared 
separation. That is, the second mode natural frequency is four times the fundamental, the 
third mode is nine times the fundamental, etc. The same structural element in stretching 
or twisting tends to have an integer spacing such as 1, 3, 5, ... for the rod fixed at one end 
and free at the other. As an example of the engineering approximation for a damped beam 
in bending, consider the simply supported uniform beam with a damper of rate c at x — £. 
For such a beam the frequencies and mode shapes are given by 


u>n — 


= 


/n7r\2 El 




n = 1,2, . 


w n {x) = sin x 


and, as is always the case, 


2C ri u n ~ 


■^^Effective 


where 



TO 

~2 
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. du dv dw 

A = h 1 

dx dy dz 

Introducing the displacement vector q and the unit vectors i, j, and k, 


and noting that 


and 


q = u{x , y, z)i + v(x, y, z)j + w(x, y , z)k, 


A = divg 


— t A dA * dA-* dA r 

gradA “ 8? ' + V + di k ' 


A = 


Ev 


G 


E 


(1 + i/)(l — 2i/)’ ~ 2(1 + 1/)’ 

the foregoing scalar Navier equations are combined into the single vector equation 


(A + G)gradA + GV z q — pq + aq. 
In free, undamped oscillations when a = 0, q — qe^ nt , and 




(A + G)gradA n + GV 2 q n = -pu 2 q n 


,2s 


Now taking the damped oscillation solution in the form 


q = q n z Xn \ 


there results the characteristic equation 

+ ('£») + u 2 = 0 . 

It follows then that the damping fraction for the n th mode of oscillation is given by 

«•-(=)©■ 

In the cases of spot damping or concentrated and localized damping, it is intuitively obvious 
that the effective mode damping fraction will exceed the result above, provided that 
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the total distributed pervasive damping equals the concentrated damping in magnitude, and 
that an anti-nodal region is employed for the damping treatment. Needless to say, dual 
modal degeneracy is excluded in this observation, where it has been seen that one of the two 
degenerate modes can be undamped in such cases. 
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6.3 Examples of Damped Beam Vibration 


The governing differential equation for a viscously damped beam in transverse bending 
vibrations is 


EI 


d^w 


( d 2 w dw\ 
dx 4 ~ ~ + a ~dt ) * 


The beam is assumed to be uniform in its flexural rigidity El and mass per unit length 
fx. The damping constant per unit length cr will take several forms in the examples which 
follow, ranging from a constant to a single or multiple discrete dampers. In every case 
an energy relationship which is now deduced provides the exact solution for the damping 
fraction when the exact, damped eigenvectors are known. It is also the basis for an acceptably 
accurate engineering approximation when the damped eigenvectors are approximated by the 
undamped ones. 

Taking the exact solution in the form 

w(x, t) = w n (x)e^ nt , 


the partial differential equation is reduced to an ordinary differential equation which follows 
below as 


El 


, v2 


dx 4 


+ (j,X w n + aX n w n = 0, n — 1,2, 


Since these eigenvectors w n are, in general, complex vectors, multiply the foregoing equation 
by the complex conjugate eigenvector and take the definite integrals over the beam span 
from x = 0 to x = l. This results in the equation 

rt dfiw n tl d 

EI Jo IbF™*’ dX + ***“ Jo + An Jo dx = °* 

Integrating by parts and invocating the beam’s boundary conditions yields the algebraic 
equation 


+ 


d r-\ f EI Jq 


fo <jw n w* n dx 


fxfo w n w* dx J 


VJL mm 


Xn + 
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Defining 


•M?) 


'Jo&jfe 

m/o ™nWn dx 


then 


■-(m s 


Jp a WnWn dx' 
w n Wn dx _ 


In the case of pervasive damping when a is a constant, 

Cnpervasive = (^) (“) > 71 = 1,2,.... 

More generally, when a is not constant over the span and the modes have the usual numerical 
spacing, w n is approximated by w nu , the undamped modal vector. Accordingly, 

,2 


„2 a, ( Jo (llH 

n VWl 


and 


C ” ~ (2 u n 

For example, consider the case of the si 
this case the exact solution follows as 


Jo ° w lu dx 

dx, 

supported beam with pervasive damping. In 


in 


\2u n j \nj \ i ) y 11 


7b 1, 2, ... , 


so that 


Cn — 


iC 


:, n = 1,2,... . 


2nn^y/JbEI * 

Consider the case now of three discrete dampers of rate (C/3) located at x = 1/ 4, i/2, and 
3^/4. Approximating w\ by the undamped mode 

Cl 


7 TX 


wi a = sin— , , Cl = 


3ttVa®7 


(i) +si “ 2 (f) +s “ 2 (t)] ’ 


and 


, c (1_\ _9L 

U \37T 2 / y/JIET 
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This approximation for the damping fraction is 


Ci “ (— ) 

\2viJ [m/o ™i u dx - 

The fundamental mode of the cantilever is now approximated by a polynomial approximation 
which satisfies the undamped cantilever beam boundary conditions. This is 


Computing the various terms in the approximation and comparing the results to the beam 
tip mounted damper yields the tabular comparison below. 



6.4 Bending Vibration of a Discretely Modelled Damped Beam 

Consider a cantilever beam of length i and constant flexural rigidity El. Concentrated 
masses mi, m 2 , and m 3 are attached at x = C/3, x — 2C/3, and x = t as illustrated in 
Figure 9. A damper of rate G is placed at x = l. Taking mi = m 2 = m 3 == m and neglecting 
the distributed mass of the beam compared to the effects of the concentrated masses, the 
three coupled equations of motion in matrix format are 

1 0 0] (toil fc \ [0 0 0] wi 

0 1 0 < *2 > + ( — ) 0 0 0 <W2 

0 0 lj [ibz ) [0 0 lj (ws 



51 


PRECEDING PAGE BLANK NOT FILMED 




A detailed numerical comparison is now carried out between the fundamental mode 
damping fraction approximation and an exact solution employing a digital computer code 
for the damped eigenvectors and associated complex eigenvalues. These digital computer 
results are presented in Appendix C. 

The data employed is as follows: 


mi = m2 = m3 = m= 1.0, 


and 


166.154 -95.538 24.923 

[^Effective] = —95.538 91.385 —33.231 

[ 24.923 -33.231 14.538 

the damper constant C is varied to determine the influence of the magnitude of C on the 

accuracy of the approximation with 

C = 0.00, 0.10, 0.20, 0.30, 0.40, 0.50, 1.00, 1.50, 2.00. 


The engineering approximation employs a Dunkerley type approximation followed by a 
matrix interaction to approximate the fundamental mode natural frequency and undamped 
modal pattern employing a pocket type calculator. This approximation yields the following 
data: 

u>i = 0.8775 (radians/second), 


r n>i 1 


f 0.1565 ' 

< w 2 . 

► ^ < 

0.5318 ► 

>3 J 


[ 1.000 y 


and the formula 


Ci = 0 . 436 a 


This permits a direct tabular comparison which follows below and which illustrates the 
excellent accuracy afforded by the method of approximation. 
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c 

“Approx.” 

“Exact” 

% Error* 

0.00 

0 

0 

0 

0.10 

0.0436 

0.043608 

0.018 

0.20 

0.0872 

0.087227 

0.031 

0.30 

0.1308 

0.130851 

0.0390 

0.40 

0.1744 

0.174499 

0.0568 

0.50 

0.2180 

0.2181739 

0.0798 

1.00 

0.4360 

0.437176 

0.2697 

1.50 

0.6540 

0.657797 

0.5806 

2.00 

0.8720 

0.8807479 

1.003 


* Percent error is defined as “E” 



It is seen that for practical engineering levels of damping, the percent error is a small 


fraction of one percent. Even at levels approaching critical damping the error is still only of 


the order of one percent. 


6.5 Coupled Bending and Torsion Vibrations of a Discretely Modelled 
Damped Beam 


Now consider the system illustrated in Figure 10. This six degree of freedom system 
extends the previous illustration by introducing coupled bending and torsional oscillations. 
The rationale is that torsional modes are, generally speaking, not as widely separated 
numerically as bending modes, so that the method of approximation might yield larger 
errors as the damper constant is increased. Once again 


m\ = m 2 = m 3 = m = 1 , 


but mass moment of inertia effects are introduced with 


h = J 2 = J 3 = 1 . 


The offset of m 3 is taken to be = 1 . Also 


(GJ/L)i = ( GJ/L) 2 = { GJ/L) 3 = 1 


53 




and 


(El/tf) i = ( EI/I?) 2 = ( EI/L% = 9. 

To test the hypothesis that the engineering approximation will still yield results of 
acceptable engineering accuracy, the third of the six modes of damped coupled bending 
and torsion is examined. The natural frequency and associated undamped modal vector are 
computed and are as follows: 


u>z — 1.341 (radians/second) 


fw 1 > 


' 0.1683 ) 

h 


0.7863 

I W2 

] h 

> S i 

0.5504 
0.1592 ' 

| 


1.0000 

l #3 j 

3 

, -.7540, 


The approximation for £3 is found to be C3 “Approx,” ~ -2249(7. The exact digital computer 
results are presented in Appendix D and yield the following tabular comparison: 


C 

^3 “Approx.” 

£3 “Exact" 

% Error 

0.10 

0.0225 

0.0233 

3.56% 

0.20 

0.0449 

0.0469 

4.45% 

0.30 

0.0675 

0.0709 

5.03% 

0.50 

0.1124 

0.1213 

7.92% 


Here the percentage errors are significantly larger than for bending vibrations only, but 
are still of acceptable engineering accuracy for preliminary design purposes when the type, 
location, and quality of damping are of principal interest in the system design. 

6.6 Vibration of a Plate With Spot Damping and Different Modes 
With Matching or Nearly Matching Natural Frequencies 

Figure 11 shows a simply supported uniform rectangular plate with a concentrated 
damper force at the interior point (x,y). The governing differential equation of free motion 
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IS 


Z)V 4 u/(a;, y, t ) + pw(x, y, t) + dw(x, y, t ) = 0. 


The undamped modal patterns are readily seen to be 

, x /rmtx\ . /"rmy\ 

w mn\ x i y) = sm — J sm (-j— J . 

The associated natural frequencies are 

"L = QV + JvV) 2 (£), 

where N = (|) . Now consider the case where differing modal patterns w rs have the same 
natural frequency. The plate aspect ratio N is then related to the modal integers m, n, r, 
and s by 

N= |(r 2 -m 2 )/(n 2 -s 2 )| 1/2 . 

It is shown analytically in References 2 and 3 that this “spot” damping treatment fails 
due to the degeneracy of modes mn and rs when their frequencies match. In effect a nodal 
point of the damped mode occurs at the damper location (x,y). 

A numerical approach also provides a demonstration of this anomaly. A modal expansion 
of undamped plate modes is employed followed by a numerical solution of the characteristic 
determinant and polynomial by computing the damped characteristic values as the plate 
aspect ratio is systematically varied to produce two modes whose natural frequencies are 
close to one another and ultimately match. The first nine modes are coupled by the damper 
as the aspect ratio is varied for several nominal damping levels ranging from 5 percent to 50 
percent of critical in the fundamental plate mode. The frequency match occurs in the fourth 
and fifth modes. The table which follows below is for a plate of constant area of one square 
meter as the aspect ratio is varied. The damper is placed at (x/a) = ( y/b ) = (l/7r). It is seen 
that the loss of damping is almost complete over a significant range of frequencies above and 
below the match at aspect ratio N = .7744 and frequency ratio p = 1 (o = .880, b = 1.1364). 
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Table: A — ab = 1.0 Meters Squared; (x/a) = ( y/b ) = (I/ 71 -) 


a 

P 

N = a/b = o 2 

lO 

o 

II 

< 4-1 

£ 

*ref = -10 

5 re f = -25 

^ref — 

.65 

0.5529 

.4225 ^ 

L70U5 

3.3506 1 

8.7111 

12.2702 

.70 

0.6372 

.4900 

1.7019 

3.3618 

8.7960 

10.2591 

.75 

0.7307 

.5625 

1.7029 

3.3696 

8.9389 

0.2759 

.80 

0.8316 

.6400 

0.0423 

0.0834 

0.1856 

0.2574 

.83 

0.8937 

.6889 

0.0387 

0.0743 

0.1392 

0.1346 

.85 

0.9345 

.7225 

0.0357 

0.0616 

0.0726 

0.0379 

.86 

0.9572 

.7396 

0.0325 

0.0469 

0.0328 

0.0071 

.87 

0.9785 

.7569 

0.0224 

0.0193 

0.0045 

0.0052 

.88 

0.9998 

.7744 

0.0000 

0.0000 

0.0001 

0.0002 

.89 

1.0211 

.7921 

0.0205 

0.0196 

0.0160 

0.0170 

.90 

1.0423 

.8100 

0.0287 

0.0401 

0.0404 

0.0406 

.91 

1.0635 

.8281 

0.0306 

0.0500 

0.0632 

0.0658 

.93 

1.1058 

.8649 

0.0309 

0.0564 

0.0942 

0.1170 

.95 

1.1478 

.9025 

0.0302 

0.0572 

0.1097 

0.1437 
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7.0 A SMART DYNAMIC VIBRATION ABSORBER: A 
COLLATERAL DAMPING APPLICATIONS TECHNOLOGY 


The dynamic vibration absorber is a subsidiary dynamic system to be attached to 
the primary system. Typically the primary system is exhibiting an undesirable dynamic 
response to a simple harmonic forced vibratory excitation, usually inherent in the system 
and not subject to significant frequency or magnitude changes to reduce the undesirable 
response. Accordingly the absorber is tuned to the offending forcing frequency and generates 
an opposing force, thereby reducing the response to zero (i.e., enforcing a node) or to an 
acceptable level. However in the absence of damping in the absorber, two new responses 
result at neighboring frequencies near the original offending one. When damping is 
introduced into the absorber these neighboring responses can be significantly attenuated. 
However this results in a significant reduction in the efficiency of the absorber compared to 
an undamped one as illustrated in Figures 12 and 13 (cf. References 4, 5, and 6). 

A smart dynamic absorber would be one which benefits from damping at the “side” 
frequencies by attenuating these new responses to a negligible level, while ignoring the 
presence of damping at its primary frequency, thereby having the potential to have the 
“best of both worlds” : full vibration suppression at the primary frequency with negligible, 
damped response at the “side” frequencies. This can be accomplished via the principle of the 
“notch filter.” In effect the damper valve is sharply frequency sensitive so that it produces its 
normal, large damping magnitude, except over a very narrow range of frequencies centered 
at the primary excitation frequency. A classical dynamic system where this is the case is 
the “Bridge-T Network” familiar to electrical engineers. Figures 14 and 15 illustrate such 
a network and its frequency response characteristics. Figures 16, 17, and 18 illustrate a 
purely mechanical dynamic system with similar characteristics. It should be recognized that 
the frequency sensitive damper valve results in a closed loop dynamic system. As in all 
such systems an engineering “trade-off” must be examined. The “loop gain” and optimum 
efficiency of the “smart” dynamic absorber as a closed loop system must be weighed against 
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the transient response characteristics of the new system. Generally the transients can be 
expected to be more “skittish” with the “smart” absorber. In fact, dynamic instability can 
ensue in the closed loop system if insufficient care is taken in adjusting the loop gain of the 
valve-absorber-primary system. 
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8.0 DYNAMIC STABILITY 
BOUNDARIES FOR BINARY SYSTEMS 

The notion of a quartic damping interaction polynomial equation developed above can 
be extended to include the dynamic stability boundaries for binary systems (Ref. 7 ). That 
is systems which can be characterized by two degrees of freedom and in which both energy 
input and energy dissipation are factors. In such systems, the various parameters, especially 
the frequency proximity of the pair of oscillators comprising the system, determine the 
boundaries of dynamic stability and the “trade-offs” useful in preliminary design analysis. 

Dynamic stability boundaries are now developed for linear, two-degree-of-freedom sys- 
tems with damping and elastic couplings. Special emphasis is placed on the influence of 
natural frequency proximity and those instabilities which stem from skew symmetric stiff- 
ness properties. These arise in numerous engineering and physical systems, but especially 
in aeroelasticity and flight dynamics, as in the case of wing flutter and aircraft stability and 
control characteristics, respectively. New insight is provided into the destabilizing effects of 
the dreaded modal “resonance.” 

The generalized co-ordinates of the binary system are denoted by q\ and qi- These 
represent modes of undamped free vibration, the natural frequencies of which are ui\ and u> 2 > 
and the generalized masses of which are and M22, respectively. The damping matrix 
is symmetric and positive definite with elements C\\, C\2, and Cfo- The stiffness matrix is 
skew-symmetric and positive definite with elements Ku, K\2, and K22, and K\\ = Mnuif, 
K22 = M22&2- The equation of the system is 

Mil 0 fgil On C \ 2 ( q\ 1 Ku K t2 f 9i 1 _ Q 

0 M 22 [ 92 / _Ci2 C 22 \<12) ~Ki2 K 2 2 l^J 

The system characteristic values A can be developed as follows. Let 
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and refer these characteristic values to the natural frequency uq, where rj = (A/oq). 
Upon introducing the dimensionless parameters p = (C 22 /Cii)/(M 2 2 /Mu) > 0, ( = 
(Cu/2M n un) > 0, S 2 = (C n Ci2 - Cf 2 ) /M n M 22 u; 2 > 0, T 2 = K 2 2 /M u M l2 uf > 0, 
and r =ui 2 /u 1, the system characteristic equation becomes 

77 4 + 2£(1 + p)rp + (1 + r 2 + 8 2 )rj 2 + 2 £(r 2 + + (r 2 4- T 2 ) == 0. 


At neutral stability a sustained oscillation will occur at dimensionless frequency v, so that 
77 = It follows at once that 

v 2 = (r 2 + p)/(l + p) and i/ 4 — (1 + r 2 + 8 2 )u 2 + (r 2 + T 2 ) = 0. 

Solving for T 2 as a function of S 2 , r 2 , and p gives 




/l+r 2 + «5 2 A 

V 1 + P ) 


Upon introducing scale changes through the new variable E, where 


E = (r z - 1)/(1 + p) and r 2 = 1 + (1 4- p)E, 


and noting that 

i/ 2 = l + £, 

the stability boundaries take the simple parabolic form 

(T 3 /p) = E 2 + (6 2 /p)(l + E). 

In this form the destabihzing and stabilizing parameters T 2 and 8 2 are related via the 
new frequency proximity variable E. The dreaded “modal resonance” phenomenon is now 
especially transparent. When the natural frequencies uq and u 2 are equal, E = 0 and T = 8: 
that is, the damping characterized by 8 is needed to provide neutral stability in the presence 
of the destabilizing effects of F. Anything less results in a divergent oscillation. It is also 
seen that when T = 8 and E = 0, the frequency of the sustained oscillations is uq = u 2 . 
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More generally, the frequency of the sustained oscillations at neutral stability differs from 
u>l when r/1. In Figures 19 and 20 are shown the stability boundaries and instability or 
“flutter” frequency, respectively, as functions of the frequency proximity variable E. 

The physical significance of the results illustrated in Figures 19 and 20 can be considered 
by an example calculation. Suppose the modal damping parameters C\\ and C 22 and the as- 
sociated generalized mass parameters Mu and M 22 are such that p = (C 22 /C 11 ) / (M 22 /M 11 ) 
= 0-50. Also, the consolidated damping-energy dissipation parameter S 2 = 0-20 and the 
modal frequency ratio parameter U 2 / 0 J 1 = 1-10. The frequency proximity variable E then 
is 0-1400. It follows at once that V 2 / p = 0-4756 and, for p = 0-50, T 2 = 0-2378. Accord- 
ingly, neutral dynamic stability results when F 2 = K 2 2 /MhM22^>i = 0-2378. When K \2 > 
0-4876(MnM22) 1 / 2 wf, a divergent oscillation occurs. The frequency of this instabihty, the 
“flutter” frequency, follows as illustrated in Figure 20 as a;g utter = l-068u;i. Since also 

r 2 = Kl 2 /M 11 M 2 24 = (^ 2 /^ 11 ^ 22 ) (u; 2 M ) 2 = 0-2378 

and (^/^l) 2 1*21, it follows that dynamic instabihty occurs for K\ 2 j (K 11 K 22) 1 / 2 > 0.4433. 

In conclusion, it is to be noted that the binary system stability boundaries are universal 
ones. They yield the levels of the destabilizing parameters in terms of generalized damping 
and frequency proximity parameters. It is also clear that the modal natural frequency ratio 
is of crucial importance in binary system dynamic stability. Severe and perhaps unstainable 
damping requirements result when the two natural frequencies match or nearly match. Hence 
the concept of the “dreaded modal resonance.” 
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9.0 CONCLUSIONS 


Methods for approximating the effects of viscous or equivalent- viscous type damping 
treatments on the free and forced vibration of lightly damped aircraft-type structures have 
been developed. Similar methods have been developed for dynamic hysteresis- viscoelastic- 
type damping treatments. In all cases it is clear that these relatively simple, energy- 
based methods yield acceptably accurate engineering approximations for preliminary design 
purposes and in most instances can supplant much more complex, finite element-type digital 
computer computations. 

Selected illustrative computational examples for a variety of structural elements have been 
carried out. This type of computation should be continued and extended to illustrate the 
procedures and methodology for an entire resonating airframe. This should be accompanied 
by suitable experimental validation as well as comparison with the digital computer, finite 
element approach. 

It is noteworthy that in the case of lightly damped structures, the apparent complication 
of intermodal coupling due to damping can be neglected with one exception. In the case 
of differing modes having matching or nearly matching natural frequencies it is necessary 
to carry out a relatively simple interaction calculation which determines the distribution of 
damping between such modes. 
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Figure 1.- Equivalent viscous damping as a function of loss factor 



Figure 2. Elastic rod embedded in viscoelastic material. 
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Figure 9.- Discrete model of a cantilever beam with a viscous damper 
at the free end. 
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Figure 10. - A discrete model of coupled bending and torsion with a viscous 
damper for bending displacement. 
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Figure 11.- Rectangular plate with spot damping. 
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Figure 16.“ Mechanical analog of notch, filter smart dynamic absorber. 
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Figure 17.- Frequency response of mechanical notch filter as a function of 
frequency ratio, viscous damping ratio, and feedback gain. 
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Figure 19 •- Dynamic stability boundaries for binary systems 
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Figure 20*- Binary flutter frequency as a function of modal frequency ratio. 



APPENDIX A 

A CLOSED FORM SOLUTION TO THE DAMPING 
INTERACTION QUARTIC EQUATION 


The damping interaction quartic equation can be written in the form 


f(rj) = J7 4 + ar) 3 + brj 2 + crj + d = 0, 


where 


V — > a — ^ (Cm* + C/3/?)i b — 1 + "I" 4 (CaaC/3/3 — Ca/jC/ta) 


c = 2 


[<-©’*<*] -ey 


and 


{X}-& 


,At 


The quartic polynomial may be factored into the product of two quadratic polynomials 
f(v) = fl(v)f2(v), where 


fi(v) = y 2 + 


| “ ^/(|) + 2 ^ -6 *7 + ^2/ + y/y^dj = 0 


and 


fiiv) = v 2 + 


| + y (i) +2y-b r] + {y- \]y 2 -d^ = 0. 


The parameter y is known as the resolvent and is a real root of the cubic equation 

F(y) = 8y 3 — 46y 2 + 2(ac — 4d)y - ^c 2 + d(a 2 — 4&) j = 0. 

This approach is due to Ferrari. The cubic itself also has a closed form solution due to 
Tartaglia and Cardan, thereby completing the closed form solution to the quartic equation. 
Rewrite the cubic equation for the resolvent in the form 


F(y) = y 3 + By 2 + Cy + D = 0, 
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where 


B = -h, C rn J(ac - 4 d), D = -i [c 2 + d(a 2 - 46)] . 

The cubic equation can be reduced to one with the second degree term absent by the change 
of variable x ■— y — ^ . This results in the reduced cubic 

F(x ) = x 3 4- px + q = 0, 


where 


_ B 2 , _ SC 2B 3 

p = c ~ T and , = 


Changing variables once again, let x = z — This yields a quadratic equation in the variable 


(z 3 ) 2 + q(z 3 )- |^=0. 


Solution of this quadratic equation yields 


— f±vS, 


where 


p 3 g^ 

fls fr + T' 


The three cube roots of 2 3 yield x±, x% and x%: 


In turn the three roots of the resolvent cubic follow: 

p b p b p b 

w = 21 “ 3^ “ 3’ m = Z2 ~^~y M = Z3 “3^“3- 

In cases of practical engineering interest y\ will be the real root of interest in determining 
the effective damping in each of the interacting modes whose natural frequencies are close 
to one another. Then from /i(p) and above and neglecting the very small term 
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4(CaaC/3/3 ~ C apCpa) compared to unity* 


1 

2 


In the case of a frequency match when 


f 




< (Caa + C/3/3) ± ^ 

(Caa + C/3/3) 2 + %1 ~ 

i+i^i 

\UaJ 

> 

4 


J 2 =1, in = i, 


and 


CEffective ~ ® 01 (Caa + C/3fj) • 

In the case of widely separated natural frequencies where | < 1 and the system is either 

spot damped or lightly damped 

CEffective - \ [(Caa + C/3/3) ± \J (Caa + C/3/3 ) 2 + ( 2 Vl ~ !)] • 

Since neighboring modes will then be effectively uncoupled by damping 

CEffective = Caa or C/3/3- 

This implies that 


(Caa + C/3/3) 2 + (2y — 1) = (Caa “ C/3/3 )“* 


Solving for yj when 


® <l ' !' = 5 (1 - 4C ““ < » ) - 

Thus it is seen that the resolvent varies over the range 


«(! “ 4 CaaC/3/3) ^ 2/1 ^ 1» 


as varies over the range 


o<f ^) 2 < 1 . 

\u a J ~ 


‘This is exactly zero in the case of spot damping. 
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APPENDIX B 

AN APPROXIMATIVE SOLUTION TO THE DAMPING 
INTERACTION QUARTIC EQUATION 

It has been shown in Appendix A that the effective damping fraction is given by 

(Caa + Cpp) ± 


CEffective ~~ 2 ' 


N 


(Cao+C^ + ZOT- 1+ (^) 




where y\ is the real root of the resolvent cubic equation. It has also been shown that as 
) varies from zero to unity, the extremes of widely separated and matching natural 

\U)aJ 

frequencies of neighboring modes, the resolvent varies in the range 

^0- - ^CaaCpp) < VI g 1- 

Now consider the mean damping fraction C, given by 


and the ratio r given by 


C = ~(C cm* + Cpp) 


( (Effective^ 

'"ct - ;• 


The foregoing equation for the effective damping fraction can be rewritten as follows: 


k 

1- 

l N 

L 


i A 


1 + 


yftt l yi 
4£ 2 


It is seen that as j — > 1, 2/1 — ► 1, and r —* 0 or 2. It is also seen that as 0, 

m -*■ 5 (i -KaaCpp) and r — ► or . (Effective approaches C aa or Cp/3 of practical 

engineering interest in the range of modal frequency ratios near a frequency match. In this 


case 


where 


( 


2 V- 1 -* 

u a J 


0 < 6 * < 1. 
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r — \ 


r 


(1 — yi) - (f ) 



1-2 

(Caa + Cpp) 2 

> 

J 


For 6 2 <C 1, the binomial expansion gives the approximation for r: 


L 

'i— * 
i 

*5 

1 

^5? 

1 


(1 - yi) - (§) 

A 

l 2 _ 

(1 - yi) - j 


r 

(Caa + Cpp) 2 

J 


(Caa + Cpp) 2 

1 i 

[ 

(Caa + Cpp) 2 

J 


Approximating the resolvent y\ as 


V 1 




6 2 \ 2c4 

7 O 


yields the approximation for r as 

7# 2 


r = 


2 / -o \ 2 

and 2 — 


0 


7# 2 


v Caa + Cpp) V ^ aa + ^ > 

7 may be approximated by reference to several numerical solutions to the interaction quartic. 
Referring to the data in Figures 4, 5, and 6, 


A geometric interpretation of the damping interaction between neighboring modes is also 
possible as follows: rewrite the Ferrari closed form solution for the effective damping fraction 
as _ 


(r-iy + 


introducing the new variable 


q = ± 


(Caa + Cpp) 2 


1 + (S) ~ 2yi 

(Caa + Cpp) 


i; 


yields a family of ellipses whose semi-major axis depends on the damping fraction levels of 
the neighboring modes; the general equation is given by 


(r - l) 2 + q 2 = 1; 
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this is seen to be a circle of unit radius with center at abscissa q = 0 and ordinate r — 1, 
Approximating the resolvent y\ or solving the resolvent cubic precisely yields the value of q 
and, in turn, leads to values of r. 



APPENDIX C 

DIGITAL COMPUTER SOLUTION TABULATION 
FOR DAMPED BENDING VIBRATION 
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C = Q 

SYSTEM EIGENVALUES DAMPINS 

REAL IMAGINARY PERCENT CRITICAL 


O.OOOOOOOOE+OO 0. 1 5436864E+02 0, OOOOOOOOE+OO 
O.OOOOOOOOE+OO -G.15436864E+02 0. OOOOOOOOE+OO 
0, OOOOOOOOE+OO O.57454960E+OI 0, OOOOOOOOE+OO 
0, OOOOOOOOE+OO -0. 57454960E+01 0* OOOOOOOOE+OO 
O. OOOOOOOOE+OO -0. 87721 309E+00 0. OOOOOOOOE+OO 
0, OOOOOOOOE+OO .0, 87721 309E+00 0, OOOOOOOOE+OO 


SYSTEM EIGENVECTORS 


RON COL 1 2 

REAL I HAG I NARY REAL IMAGINARY 


1 O.iOOOOOOE+Ol O.1064962E-27 OJOOOOOOE+Ol -0.46S1889E-27 

2 -0.6989843E+OO 0.2269553E-26 -0.6989843E+00 -0.54L3558E-27 

3 0.21 51 9 1 3E+00 -0 , 1 01 053 IE- 26 0.21519138+00 Q.5053640E-28 

4 0.4986863E-25 0.L543686E+02 0.9043501E-26 -0.1543686E+02 

5 -D. 2343B64E-25 -O.1079O13E+O2 -0. 1047923E-25 0. 1079013E+02 

6 0. 1076006E-25 0.3321879E+01 0.4132842E-26 -0.3321879E+01 


SYSTEM EIGENVECTORS 


RON COL 3 4 

REAL IMAGINARY REAL I HAG I NARY 


1 0.841 71 92E+00 0. 1289842E-21 G.8417192E+00 -0.1291827E-21 

2 0. lOOOOOOE+Ol “0. 46347 16E-25 O.IOOOOOOE+Ol 0.1565100E-25 

3 -0. 663293 1E+00 -0.7749615E-21 -0.6632931E+00 0.7725248E-21 

4 -0.3344544E-21 0.4836094E+01 -0.3350736E-21 -0. 4836094E+01 

5 -0.5011495E-21 0.5745496E+01 -0.5044752E-21 -0.5745496E+01 

6 0.2752724E-20 -0.3810948E+01 0.2746154E-20 Q.3810948E+01 


SYSTEM EIGENVECTORS 


RON COL 5 6 

REAL IMAGINARY REAL I HAG I NARY 


1 0.1564141E+00 -0.1596148E-25 0.1564141E+00 0.3077212E-25 

2 0.5316363E+00 -0.17323O2E-25 Q.5316363E+00 0.3347339E-25 

3 O.IOOOOOOE+Ol -0.3175165E-28 O.IOOOOOOE+OI -0.3549874E-29 

4 0 . 402 1 223E-25 -0 . 137 2085E+00 O.7748807E~25 0.1372085E+00 

5 0,21 3841 4E-24 -0.4663584E+00 0.4I21986E-24 0.4663584E+00 

6 0.4810201E-24 -0.8772131E+00 0.9268212E-24 0.877213IE+00 
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C = .1 


SYSTEM EIGENVALUES 

DAMPING 

REAL 

IMAGINARY 

PERCENT CRITICAL 

-O.15O84290E-O2 

-0. 15436854E+02 

0.977I6086E-02 

-O. L508429QE-02 

O.15436854E+02 

O.97716086E-O2 

“0. 30237098E-O3 

-0,57453482E+0i 

0. 178i8O36£+00 

-0. IO237O99E-01 

0,574534B2E+0i 

0.17818034E+00 

-0.38254473E-01 

0.87 640033E+00 

0.43G08019E+O1 

-0.38254473E-01 

-O.87640033E+00 

O.436080 3 9E+01 


SVSTEH EIGENVECTORS 


RON COL 1 2 

REAL l HAG INARY REAL I MAS I NARY 


\ 0. iOOOOOOE+Oi -0.2283601E-17 O.iOOOOOOE+Oi -O.2543099E-37 

2 ~0,6989836E+QO 0. 1077038E-03 -0.6989836E+00 -O. 1O77038E-03 

3 0. 2151 B16E+00 -0. 1455723E-02 0.2151816E+00 0. 1455723E-02 

4 -0, J508429E-02 -0., 3543685E+02 “0.1508429E-02 0. 1543685E+02 

5 0.2716975E-02 0. 1079011E+02 0.2716975E-02 -0. 10790UE+02 

6 “0.2279633E-01 ~0.332i725EKii -0, 2279638E-01 0.3321725E+01 


SYSTEM EIGENVECTORS 


RON COL 3 4 

REAL IMAGINARY REAL IMAGINARY 


1 0/8416667E+00 -0.3044757E-02 0.8416667E+O0 0.3044757E-02 

2 0. 1000000E+01 0, 1902775E-36 0. lOOOOOOE+Ol -0. 16I5461E-1G 

3 -0.6&3O879E+0O 0.122933BE-01 -0.6630879E+0Q -0.3229338E-O1 

4 -0.2630943E-01 -0.4835637E*O3 -0. 2610941E-01 0. 4S35637E+01 

5 -0. 102371 OE-Oi -0. 574534SE+01 -0.1023710E-01 0.5745348E+03 

6 O.7741785E-03 0.3809545E+01 0, 774i735E-*01 -0.3809545E+01 


SYSTEM EIGENVECTORS 


RON COL 5 b 

REAL IMAGINARY REAL IMAGINARY 


1 0.I563800E+00 0.7433139E-03 0, 156380GE+00 -0.7433139E-03 

2 0, 5315824E+00 0, 1177009E-02 0. 531 5826E+00 -0.11 77009E-02 

3 0.1000000E+03 0. 2943609E- I 6 O.IOOOOQOE+Oi -0.8267042E-I8 

4 -0.6633677E-O2 0. 1370231E+00 -0.6633677E-02 -0. I370231E+00 

5 -0.2t3£494E-0i 0.4658341E+00 -D.2136694E-01 -0.4658341E+0O 

6 -O.3825447E-01 0.8764003E+00 -0.3825447E-03 -0.8764003E+OO 
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C = .2 

■SYSTEM EIGENVALUES 


REAL 


-0.3QU4747E-02 
-0.301 6474 7E-02 
-0, 204630 16E-01 
-0. 204630 16E-01 
-0. 765205 10E-0 5 
-0. 765205 10E-01 


IMAGINARY 


0.15436825E+02 
“0. 15436825E+02 
-0.57449049E+01 
0.57449049E+01 
"0.B7395648E+OO 
0.S7395648E+00 


DAMPING 

PERCENT CRITICAL 


0.T954O772E-01 
0.19540772E-01 
0 . 3561 919 I £+00 
0.3561 91 91E+00 
0.B7222739E+01 
0.87222739E+01 


SYSTEM EIGENVECTORS 


RON COL 1 

REAL IMAGINARY 


1 D.IOOOOOOE+OI 0.33881 32E-18 

2 -0. 698981 5E+00 -0.2153812E-03 

3 0.215152 6 E +00 0.2911066E-02 

4 -O.3016475E-Q2 0. 1543682E+02 

5 0.5433262E-02 ~0. 1079005E+02 

6 -0.4558662E-01 0.3321263E+01 


SYSTEM EIGENVECTORS 


RON COL 3 

REAL IMAGINARY 


t 0.8415093E+00 -0. 6084581 E-02 

2 0 . 1 000 00 OE +01 O.3OOO530E-16 

3 -0.6624727E+00 0.2456817E-01 

4 -0.5217516E-01 -0.4834266E+01 

5 "0.2046302E-01 -0.5744905E+01 

6 0. 1546980E+00 0.3805340E+Q1 


SYSTEM EIGENVECTORS 


RON COL 5 

REAL IMAGINARY 


1 0. 1562779E+00 -0. 14B1613E-02 

2 0.531421 4E +0 0 -0.2346193E-02 

3 0. 1000000E+01 0.1160096E-16 

4 -0. 1325333E-01 -0. 1364667EK*0 

5 -0.42715UE-01 -0.4642597E+00 

6 -0.7652051E -01 -0.8739565E+00 


2 

REAL IMAGINARY 


0. 1000000E+01 -0. 1256319E-16 
-0.6989815E+00 0.2153812E-03 
0.2151 526E+00 -0.29U066E-02 
-0.3016475E-02 -0. 1543682E+02 
0.5433262E-02 0, 1079005E+02 
-0.4558662E-01 -0.3321263E+01 


4 

REAL IMAGINARY 


O.8415O93E+0O 0.6084581E-02 
0. lOOOOOOE+Ol 0.1818749E-16 
“0.6624727E+00 -0.2456817E-01 
-0.5217516E-01 0.4834266E+01 
“0.2046302E-01 0.5744905E+01 
0, 1546980E+00 “0.3805340E+01 


6 

REAL IMASINARY 


0. 1562779E+00 O.1481613E-02 
0. 5314214E+00 0.2346193E-02 
0. lOOOOOOE+Ol 0.2602085E-17 
~0. 1325333E-01 0.1364667E+00 
-0.427151 l-E-01 0.4642597E+00 
-0.765205 IE -01 0.8739565E+00 
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C = .3 

SYSTEM EIGENVALUES 
REAL IMAGINARY 


-CM 5237 54 4E ~02 
-0.45237544E-O2 
-0. 30666 592E-01 
-0.3O666592E-O1 
-0.il 480966E+00 
-0, 11480966E+00 


0.15436775E+02 

■-0.1543677.5E+02 

-O.5744167OE+01 

0.57441670E+01 

-O/869864A8E+O0 

0.86986468E+00 


DAMPING 

PERCENT CRITICAL 


O.29305O46E-O1 
0.29305046E-01 
0, 53386601 E+00 
O.53386600E+OO 
0. 13085085E+02 
0. 1 3085085E+02 


SYSTEM EIGENVECTORS 


RON COL 1 

REAL IMAGINARY 


1 0.1 OOOOOOE +01 0.4483i76E"16 

2 -O.698978OE+0O -0..3230059E-03 

3 0, 2 151041 E+09 0.4365647E-02 

4 -O.4523754E-02 0, 1543678E+02 

5 0.8148175E-O2 -0.1078997E+02 

6 “0,6836459E-O1 0.3320494E+01 


SYSTEM EIGENVECTORS 


RON COL 3 

REAL IMAGINARY 


1 0.8412475E+O0 -0.9M4561E-02 

2 0 « 1 OOOOOOE +01 -0.3377290E-16 

3 -0.66 1449 IE +00 0.36B0585E-01 

4 -0.7815375E-01 -0.4831986E+01 

5 -0.3O66659E-O1 -O.5744167E+01 

6 0.2317033E+00 0.3798345E+01 


SYSTEM EIGENVECTORS 


RON COL 5 

REAL IMAGINARY 


1 0.1561077E+O0 -0.22O9853E-O2 

2 0.531 153 1E+0O -0.3 49968 0E~ 02 

3 0. lOOOOOOE+Ol 0.5149960E-17 

4 -0.1984495E-01 -0.1355389E+00 

5 -0.6402575E-01 -O.4616295E+0O 

6 -0.1148097E+00 -0.8698647E+00 


REAL IMAGINARY 


0. lOOOOOOE+Ol -0.2 6264 80 E - 1 6 
-0.698978OE+OO 0.3230059E-03 
0,2 15 1041 E+00 -0.4365647E-02 
-0. 4523754E-02 -O.I543678E+02 
O.8148175E-02 0.1078997E+02 
-0.6836459E-01 -0.3320494E+01 


4 

REAL IMAGINARY 


0.8412475E+00 0.911456iE-O2 
0,1 OOOOOOE +01 -O.2057274E-16 
-0.661 449 1E+O0 -0.3680585E-01 
-0.7815375E-0 1 0.4831986E+01 
-O.3O66659E-01 0.5744167E+01 
O.2317O33E+0O -0. 3798345E+01 


6 

REAL IMAGINARY 


0. 1 56107 7E+00 0.2209853E-02 
0.531 1531 E+00 0.3499680E-02 
0. 1 OOOOOOE +01 0.3815O36E-17 
-O. 1984495E-01 0.1355389E+00 
-O.6402575E-O1 0.4616295E+00 
-0. M48O97E+0O 0.8698647E+00 
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C = .4 

SYSTEM EIGENVALUES DAMPING 

REAL IMAGINARY PERCENT CRITICAL 


-0. 60298B52E-G2 0. 15436706E+02 0.3904I991E-01 
-0.40298852E-02 -0.15434704E+02 0.39041991E-01 
-0.40B34495E-01 -0. 57431354E+01 0. 7I103430E+00 
-0.40834495E-01 0. 57431354E+01 0.7U03430E+00 
-0.15313342E+00 -0.84409429E+00 0.1744990IE+02 
-0. 15313342E+00 0.84409429E+00 0.17449901E+02 


SYSTEM EIGENVECTORS 


RQN COL 1 2 

REAL IMAGINARY REAL IMAGINARY 


1 0. 100000QE+01 0.4478332E-14 0.1000000E+01 0.4803014E-14 

2 -0, 4989732E+00 -0.4305517E-03 -0.4989732E+90 0.4305517E-03 
3- 0.2150343E+00 0.5819084E-02 0.2150343E+00 -0.5819084E-02 

4 -0.4029885E-02 0.1543471E+02 -0.4029B85E-02 -0. 1543471E+02 

5 0.1084103E-01 -0. 1078984E+02 0.1084103E-01 0.1078984E+02 
4 -0.9I12414E-91 0.33I9418E+01 -0.9U2414E-01 -0.3319418E+01 


SYSTEM EIGENVECTORS 


RON COL 3 4 

REAL IMAGINARY REAL IMAGINARY 


1 0.8408819E+00 -0.1212982E-01 0.8408B19E+00 0. I212982E-01 

2 0.1000000E+01 -0, 2549559E-14 0. 1000000E+01 -0.3174002E-16 

3 -0.4400194E+00 0.4898799E-01 -0.4400196E+09 -0.489B799E-01 

4 -0.1040020E+00 -0.4828803E+0I -0. 1040020E+00 0,4828803E*0i 

5 -0.4083470E-01 -0.5743135E+01 -9.4083470E-01 0.5743135E+01 
4 0.3082977E+00 0.3788581E+01 0. 3082977E+00 -0.378B581E+01 


SYSTEM EIGENVECTORS 


RON COL' 5 6 

REAL IMAGINARY REAL IMAGINARY 


1 0. 1558499E+00 

2 0.5307779E+00 

3 0. lOOOOOOE+Oi 

4 -0.2439459E-01 

5 -0.8528017E-01 
4 -0.1531334E+90 


-0.2922929E-02 0.1558499E+00 0.2922929E-02 
-0.4629492E-02 0.5307779E+00 0.4429492E-02 
D.4878910E-18 0.1000000E+01 0.2981554E-18 
-0.1342390E+00 -0.2639459E-01 0.1342390E+00 
-0.4579343E+00 -0.8528017E-01 0.4579343E+00 
-0.8640963E+00 -0.1531334E+00 0.8640963E+00 
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5 


C = . 


SYSTEH EIGENVALUES 

DAHPING 

REAL 

IHAGINARY 

PERCENT CRITICAL 

-0.7534486OE-02 

0, 15436617E+02 

0.48809172E-O1 

-0, 75344860E-02 

-0. 1 54366 17E+02 

0.48809172E-01 

-O.5O962256E-0I 

-0, 57418120E+01 

0.8B752904E+00 

-O.5O962256E-01 

O. 57418120E+OI 

0.88752903E+0O 

-0.19150326E+O0 

-0. 8566101 If +O0 

0/21817386E+02 

-0, 1915O326E+0O 

O. 8566101 1E+Q0 

0,2 18 l 7386E+02 


SYSTEH EIGENVECTORS 


ROW COL 1 2 

REAL I HAGINARY REAL I HAG INARY 


1 0, lODOOOOE+OI ”0,2 6264 80 E - 16 0.1000000E+01 O.5708324E-16 

2 “0. 69S9670E+OO -0.5379923E-03 -0.6989670Et00 0.5379923E-03 

3 0,21 49493E+00 0.7271004E-02 0.2149493E+00 -0. 7271004E-02 

4 -0.7534486E-02 0. 1543662E+02 -O.7534486E-02 -0.1543662E+O2 

5 O. 13571 14E-01 -0, 1078968E+02 0. 13571 14E-01 0. 1078968E+02 

6 -0.U38592E+00 0- 3318035E+01 -0. U38592E+00 -0.3318035E+01 


SYSTEH EIGENVECTORS 


ROW COL 3 4 

REAL IHAGINARY REAL IMAGINARY 


1 0, 8404135E+00 -0. 1512553E-01 0.8404135E+00 0. 1512553E-01 

2 0. lOOOOOOE+Ol 0.5597194E-17 0. lOOOOOOE+Oi -0. 2141299E-17 

3 -0.65B1876E+OO O.61O9639E-01 -0.6581876E+00 -Q.61O9639E-0I 

4 -O.1296773E*0O -0.4824725E+01 -0. 1296773E+00 0. 4S24725E+01 

5 -0, 5096226E-01 ”0. 5741812E+01 ~0,5096226E-O1 0.5741812E+01 

6 0 . 3843467E+00 0. 3776076E+01 0, 38434A7E+00 -0. 377607AE+01 


SYSTEH EIGENVECTORS 


RON COL 5 6 

REAL IHAGINARY REAL IHAGINARY 


1 0. 1555649E+0Q -G.3615637E-02 0.1555649E+00 0.3615637E-02 

2 0, 5302966E+00 -Q.5727494E-G2 O.53O2966E+0O 0.5727494E-02 

3 0» 1000000E+01 “0.2954451E-17 0. lOOOOOOE+Oi -0.4526544E-17 

4 -0.328883BE-O1 -0. 1325661E+00 -0.3288B3BE-Q1 0, 1325661E+00 

5 -0.1064598E+00 -0.453U06E+00 -0. 1064598E+00 0. 4531606E+00 

6 -0, 1915033E+00 -0. 8566101E+00 -0, 1915033E+00 O.85661O1E+O0 
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C = 1.0 


SYSTEM £] 

GENVALUES 

DAMPING 

REAL 

IMAGINARY 

PERCENT CRITICAL 

-0,1502I291E-Oi 

O, 154358 79E+02 

0,973 1408 iE-01 

-0.1 5021 291 E-01 

-0. 15435879E+02 

0,973i408t£”Ol 

-CU0054O50E+0O 

-0.573O9066E+O1 

O, 17540860E+01 

-0. 10054 05 OE +00 

O. 57309O66E+01 

O. 17540860E+01 

-O. 3844382 IE tOO 

-O, 790831 16E+00 

O.43717624E+02 

-0.38443821E+0O 

0.79088 11 AE+OQ 

0.43717424E+02 


SYSTEM EIGENVECTORS 


ROW COL i 2 

REAL IMAGINARY REAL IMAGINARY 


t 0, 1 0O0OO 0E+01 0, 2198220E-16 O.IOOOOOOE+Ol -0.1045577E-16 

2 -O.6989152E+0O -0. 1072706E-02 ~0. 6989152E+00 0. 1072706E-02 

3 0.2142262E+00 0. 1449462E-01 0.21422*2£+O0 -0.1449462E-01 

4 -0.15O2129E-OI 0.1543588E+02 -Q.1502129E-01 -0.1543588E+O2 

5 O.2705677E-O1 -0, 10788351+02 0, 2705677E-01 9.107B83SE+02 

6 -Or 2269552E+00 O.3306551E+O1 -0.2269552E+00 -0. 330655 iEtO.1 


SYSTEM EIGENVECTORS 


RON COL 3 4 

REAL IMAGINARY REAL IMAGINARY 


1 0. 8365826E+00 -0. 2964796E-01 O.8365826E+O0 0.2964796E-01 

2 O. lOOOOOOE+Ol -0.2481468E~16 0. 1000000E+01 -0.7535205E-17 

3 -0.643183OE+OO 0, U99i47E+00 -0, 6431830E+00 -0. U99U7E+Q0 

4 -0.254O201E+0O -0.4791396EKU -0.25402O1E+0O 0.4791396E+01 

5 -D. 1005405E+00 -0. 5730907E+0I -0, 1005405E+00 O. 5730907E+01 

6 0, 7518856E+00 0.3673964E+01 0.7518B56E+00 -0.3673966E+01 


SYSTEM EIGENVECTORS 


ROW COL 5 6 

REAL IMAGINARY REAL IMAGINARY 


1 0J53O536E+OO -0.6579O89E-02 0. 1530536E+00 O.&579O89E-02 

2 0.5263284Ef00 -0.1043489E-OI 0. 5263284E+00 0„ 1043489E-01 

3 O.IOOOOOOE+Ol -0.2591243E-16 0. 1000000E+0I 0.7909255E-16 

4 -Q.64O4294E-01 -0. U85180E+00 -D.6404294E-01 0. 118518OE+0O 

5 -*0. 2105935E+00 -O.4122516E+00 -0. 2105935E+00 Q.4122516E+0Q 

6 HU844382E+0O -O.7908812E*Q0 -0. 3844382E+00 0. 7908812E+00 
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C « 1.5 
'SYSTEM EIGENVALUES 
REAL IMAGINARY 


-0,224 13669E-01 
-G.22413669E-G1 
-O.1474OB69E+O0 
-0, 14740868E+00 
-0. 580 17 76 5£ +00 
-0.5801 7 765E+00 


0.15434658E+02 
-Q.15434658E+02 
-0.57132182E+O1 
0,571 321 82E+01 
0.66431856E+00 
-O.64431854E+O0 


DAMPING 

PERCENT CRITICAL 


0.1 4521 634E+0O 
0, I4521634E+00 
0, 25792756E+01 
0.25792756E+01 
0.65779726E+02 
0.45779726E+Q2 


SYSTEM EIGENVECTORS 


RON COL 


REAL 


1 


IMAGINARY REAL 


IMAGINARY 


1 0. 100000QE+01 

2 -0.4988294E+00 

3 0, 213031 1E+00 

4 -0.2241347E-01 

5 0.4037305E-01 

6 -0 . 3385404E+00 


G.3778445E-16 
-O.l 6009 24E-02 
O.2162442E-01 
O, 1543466E+02 
-0.1078616E+02 
0.3287577E+01 


0, lOOOOOOE+Ol 
-O.6988296E+O0 
0.21303.UE+00 
-0.2241367E-O1 
0. 4O37305E-O1 
-0.33854O4E+OO 


0.2574980E-18 
0, 16OO924E-02 
-0.21 624 42E-01 
-0,1 5434 66E+02 
0* 10786 16E+02 
-0.3287577E+01 


SYSTEM EIGENVECTORS 


RON COL 


4 


REAL IMAGINARY REAL IMAGINARY 


1 0, 8 3 04 75 9E +00 

2 0,1 OOO CODE +01 

3 -0. 6191851 E+OQ 

4 -0.3681 931 E+00 

5 -0.1474087E+00 

6 0. 1087446E+01 


-0.43O1845E-O1 
0.1897354E-16 
0, 1743629E+00 
-0.4738349E+01 
-0.5713218E+01 
0.351 1837E+01 


0.8304759E+00 
O, 1 OOOOO OE +01 
-0.6 191 85 IE +00 
-0. 3681 931 E+00 
-0.1474087E+00 
0.1087446E+01 


0.4301845E-01 
-0.8605855E-18 
-0. I743629E+00 
0.4738349E+01 
0.57 132 I8E+01 
-O.3511837E+01 


SYSTEM EIGENVECTORS 


RON COL 


6 


REAL IMAGINARY REAL IMAGINARY 


1 0. 1489885E+00 0.8088565E-02 

2 0.5198839E+00 0. 12856 18E-01 

3 0. lOOOOOOE+Ol -0.1395910E-17 

4 -0.9181317E-01 0.9428300E-01 

5 -0.3101656E+00 0.3379096E+00 

6 -0.5801776E+00 0.6643186E+00 


0.1489885E+00 -0.8088565E-02 
0.5198B39E+00 -0. 12856 18E-01 
0. lOOOOOOE+Ol 0, 168051 3E- 16 
-0. 91 B1317E-01 -0 . 9428300E-0 1 
-0.3101656E+00 -0.3379096E+00 
-0.58O1776E+OO -0.6643186E+00 
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C - 2 


SYSTEM EIGENVALUES 
REAL IMAGINARY 


-0. 29666716E-01 
-0.29666716E-01 
-0.19036224E+O0 
-0,19G36224E+0D 
“0. 77997104E+00 
-0, 77997 1O4E+00 


-0. 15432970E+02 
0.15432970E+02 
-Q, 56894674E+01 
0. 5689467 4E+01 
O.41939680E+0O 
-0.41939680E+O0 


DAMPING 

PERCENT CRITICAL 


0. 1922291 1E+00 
0. 1922291 1E+-00 
0.33439998E+01 
0.33439998E*D1 
0.88074793E+D2 
0.88074793E+02 


SYSTEM EIGENVECTORS 


RON COL 1 

REAL IMAGINARY 


1 0. 1000000E+01 Q.2726768E-16 

2 -O. 69871 llE+OO 0.2H9552E-Q2 

3 0.2113789E+00 -0,2861587E-0i 

4 -0.2966672E-01 -0. 1543297E+02 

5 0, 534394 4E-0 I 0.1078312E+02 

6 -0. 4478987E+00 -0.3261355E+D1 


SYSTEM EIGENVECTORS 


ROW COL 3 

REAL IMAGINARY 


1 0.8224784E+O0 -0.5479O12E-01 

2 0. 1000000E+01 0.4387631E-18 

3 -O.5876064E+QO 0. 2227034E+00 

4 -O.4682954E+0O -0.4669O34E+OI 

5 -O. 19O3622E+O0 -0.56B9467E+01 

6 O. 1378922E+01 0.3300773E+Q1 


2 

REAL IMAGINARY 


0.100O0O0E+01 0,54210 HE- 19 
-0.698711 1E+0O -0.2119552E-02 
0.2113789E+O0 0.2861587E-01 
-0.2966672E-01 O. 1543297E+02 
0.5343944E-01 -0.1D78312E+02 
-O.4478987E+0O 0.3261355E+01 


4 

REAL IMAGINARY 


0.8224784E+0O 0. 5479O12E-01 
0. IGQOOQOE+Ql 0.U04531E-16 
-O. 5876O64E+00 -O. 2227O34E+0O 
-0.4682954E+O0 0.4669034E+01 
-0. 1903622E+0O 0.5689467E+O1 
0.1378922E+01 -O.33O0773E+01 


SYSTEM EIGENVECTORS 


ROW COL 

REAL 


5 

IMAGINARY 


O.6578427E-02 

0,iO48759E-Ol 

-0.3930233E-17 

0.5507235E-01 

0. 20622 15E+ 00 
0.4193968E+O0 


REAL 


0.1435474E+O0 

O.511214OE+00 

0.100O000E+01 

-O.1147218E+00 

-0.4031306E+00 

-0.7799710E+O0 


6 

IMAGINARY 


-0.6578427E-02 
~0, 1048759EO1 
0.2347298E-16 
-0,55O7235E-01 

-0. 2D622 15E+0O 
-0.4193968E+0O 


1 0, 1435474E+O0 

2 0.511 21 4OE+0O 

3 0.10O0000E+03 

4 -0.H47218E+00 

5 -0.4O31306E+00 

6 -0,77997 iOE+00 
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APPENDIX D 

DIGITAL COMPUTER SOLUTION TABULATION FOR 
DAMPED, COUPLED BENDING-TORSION VIBRATION 
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C = 0.00 


SYSTEM EIGENVALUES 
REAL IMAGINARY 


O.OO.OOODOOE+O0 
O.OO00000OE+0O 
O, 0000000 OE +00 
O.0O00000OE+0O 
Q.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 
0.O000OO0OE+0O 
O.OOOOOOOOE+OO 
0.0O000000E+0O 
O.OOOOOOOOE+OO 
0.00000 DO 0E +00 
0 < 0 000000 OE+00 


-D. 15666482E+02 
0.15646482E+02 
O. 63089979E+01 
-O.63089979E+01 
0. 18001 IB8E+01 
■-0. 18001 188E+01 
-0.1 34071 44E+0I 
0. 13407 144E+01 
-0.93963702E+0O 
0.93963702E+00 
0 . 34696696E+0O 
-O. 34696696E+00 


DAMPING 

PERCENT CRITICAL 


Q.OOOOOOOOE+OO 
O.0O0000OOE+0O 
0.0OV0O000E+00 
O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 
0* 000O0000E+00 
O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 
0. OOOOOOOOE+OO 
O.OOOOOOOOE+OO 
O.OOOOOOOOE+OO 


SYSTEM EIGENVECTORS 


RON COL I 


2 


REAL 


IMAGINARY REAL IMAGINARY 


1 0 . 1 OOO OOOE +0 1 0, 24563 96E- 19 

2 “O.38O9764E-05 0- 1 452009E-23 

3 ~0* 7 1581O9E+O0 0. 17951B0E-18 

4 0, 9274438E-03 0. 2260271E-21 

5 0.450620 1 E+00 -0.3959879E-19 

6 -O. 22577 19E+00 0. 2538452E-19 

7 “O. 1856696E-17 -O. 1566648E+02 

8 0.12H044E-22 0.596856OE-04 

9 -0.15O6O25E-17 0.1121424E+02 
ID 0.5558654E-21 -0. 1452978E~CM 

11 -O. 1 390828E-I7 -0. 7059632E+01 

12 -0.2655448E-18 0.3537O51E+O1 


0. lOOOOOOE+Ol -0.4 086 93 4 E - 1 9 
-0.38O9764E-05 -0.5987366E-23 
“0.71581O9E+00 -0.45 18921 £“18 
0.9274438E-03 Q. 153Q077E-20 
0.4506201E+OD O.2038173E-18 
-0. 22577 19E +00 -0. 1387016E-18 
-0.146740OE-16 O, 1566648E+02 
0. 1105449E-21 -0.5968560E-04 
0.1684240E-16 -0.1121424E+02 
- 0 . I 7 07 963E ~ 19 O. 1452978E-01 
-0.8570279E-17 O.7O59632E+01 
0.4852652E-I7 -O. 3537051E+01 


SYSTEM EIGENVECTORS 


RON COL 3 


REAL 


1 -0. 7238354E+00 

2 -O. 3546887E-03 

3 - O , 69 64 22 5E +00 

4 O.134O846E-01 

5 O. 10O000OE+01 

6 “0. 50653 13E+00 

7 “0.3369603E-i9 

8 - 0 . 26 8 1 87 5E - 24 

9 “O.3105567E-19 

10 0.5637236E-21 

11 0 . 23 45222E - 1 9 

12 -0. 7649766E-20 


IMAGINARY 


O.4362551E-20 
“0.1558234E-23 
0.4663644E-20 
Q.307Q908E-22 
0.53601 30 E-2 1 
0.6294844E-21 
“0.4566676E+01 
-0. 2237730E-02 
-0.4393728E+01 
0 « 8459393E-0 1 
0.6308998E+01 
-0.3195 705 E+ 0 1 


REAL 


-0.7238354E+0O 
-0.3546887E-03 
-O.6964225E+00 
0.1340846E-01 
0.1 000 0 DOE +01 
-0.5065313E+00 
0.2767746E-20 
0.3699695E-24 
0.2081 186E-20 
-0.41 82432E -22 
-O.6882143E-20 
0.207126OE-2O 


4 

IMAGINARY 


0,1379324E"21 
0.8613099E-25 
0. 6896618E-21 
0.7496324E-24 
O,50O4443E“22 
O.2163077E-21 
0.4566676E+D1 
O. 2237730E-02 
0.4393728E+01 
-0.8459393E-01 
-O.630899SE+01 
0.31957O5E+O1 
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C = 0.00 (continued) 
ROW COL 5 

REAL IHA6INARY REAL 


1 0.7960482E-01 -0. 277701 4E-22 0.7960482E-01 

2 -0.8061736E+00 -0, 23841 24E-21 -0.8061736E+00 

3 0.2468506E+0Q -0.9222276E-22 0.2468506E+00 

4 0. 1000000E+01 0.14 13639 E-25 0. lOOOOOOE+Gl 

5 0.4258956E+QQ -0.1697927E-21 0, 4258956E+00 

6 -0.43 4254 IE +00 0, 4665079E-22 -0. 434254 iE+OO 

7 0.36O5B47E-22 0.1 43298 iE+OO 0.1B97984E-24 

8 0.3462502E-21 -0. 1451208E+01 0.1731329E-23 

9 0. H7901 IE-21 0.4443603E+00 0.6143775E-24 

10 -0.2742540E-22 0.1800U9E+01 -0. U69281E-24 

11 0.2189405E-21 0.7666627E+00 0.5091998E-23 

12 -0.9348757E-22 -0.7817090E+00 -0.4768633E-24 


SYSTEM EIGENVECTORS 


ROW COL 7 

REAL IMAGINARY REAL 


1 D. 1683205E+00 O.U25507E-26 0. 1683205E+00 

2 0.7862634E+00 0. 1435964E-25 0.7862634E+00 

3 0. 5504385E+00 0.1097700E-26 0.5504385E+00 

4 0.1592064E+00 0.2585699E-25 0.1592064E+00 

5 0. 1000000E+01 0.5656 l 33E-27 0. 1000000E+01 

6 -0.754G265E+00 0.7494179E-28 -0.7540265E+00 

7 0. 1999565E-26 -0. 2256697E+00 0.B926876E-22 

8 0. 1 170591E-25 -0. 1054155E+01 0.4865417E-21 

9 0.4620358E-26 -0.73798D8E+00 0.2465933E-21 

10 0.3189522E-25 -0.2134504E+00 0. 1567947E-20 

11 0. 123251 6E-25 -0. 13407 14E+01 0.3921293E-21 

12 -O. 83556 15E-26 0.1010934E+01 -0.3798375E-21 


SYSTEM EIGENVECTORS 


ROW COL 9 

REAL IMAGINARY REAL 


1 0. 357681 9E+0Q 0.3317002E-24 0, 1576819E+00 

2 -0. 67557 55E +00 -0.2535076E-20 -0.6755755E+00 

3 0. 5336430E+00 0.I9B2123E-24 0.5336430E+00 

4 -0. 7546734E+00 -0.3080662E-20 -0.7546734E+00 

5 O.iOOOOOOE+Ol 0.U65242E-24 0.1000000E+01 

6 -0.167456BE+00 -O.45G3340E-2O -0. 1674568E+00 

7 0.1108697E-24 -0.1481637E+00 0.5520764E-25 
S -0.2230077E-20 0.6347957E+00 -0.3998230E-22 

9 Q.4162479E-22 -0.5014307E+00 0.6825982E-24 

10 -0. 1735B40E-20 0.7091 191E+00 -0.31 222 13E-22 

11 0.8568832E-22 -0.9396370E+00 0. 1715425E-23 

12 -0.216941 IE-20 G, 1573486E+00 -0.3901747E-22 


6 

IMAGINARY 


0. 1461026E-24 
0.1199700E-23 
0.4812209E-24 
0.9138165E-26 
0.8750992E-24 
-0.2358126E-24 
“0.1 43298 IE+OO 
0.1451208E+01 
-0.4443603E+00 
-0. 18001 19E+01 
-0.7666627E+00 
0.781 70 90E+00 


8 

IMAGINARY 


-0. 132391 BE- 22 
-0.6347904E-21 
-0.2090 l 25E-22 
-0. 1302085E-20 
-0.9855082E-25 
-0. 1113612E-21 
0.2256697E+00 
0. 1054155E+01 
0.7379808E+00 
0.2134504E+00 
0.1340714E+01 
-0. 1010934E+01 


10 

IMAGINARY 


-0. 40724 16E- 25 
0.4547383E-22 
-0.937924 IE-25 
0.5538442E-22 
0.8305128E-26 
0.8085720E-22 
0. 1481637E+00 
-0.6347957E+00 
0.5014307E+00 
-0.7091 191E+00 
0.9396370E+00 
-0. 1573486E+00 
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C = 0.00 (concluded) 


ROW COL 11 12 

REAL IMAGINARY REAL IMAGINARY 


1 0.2129991E-01 0. 9425783E-20 0.2129991E-01 -0.8705053E-20 

2 0.3947948E+00 0.3412594E-18 0.394794BE+00 -0.315I543E-18 

3 0. 7420595E-01 0. 7824545E-29 0. 7420595E-01 -0.7228171E-20 

4 0.7420454E+00 0.4384725E-18 0.7420656E+00 -0.4049437E-18 

5 0, 1425404E+00 -0.U48952E-19 0, 1425404E+00 0.1041083E-19 

4 0. I OOOOOOE+G1 -0.4175283E-25 0, 1009009E+01 -0.1048417E-24 

7 0.7722742E-20 0, 7390343E-02 0.7132210E-20 -0.7390343E-02 

8 -0.3257859E-19 9.1349814E+00 -0.3008734E-19 -0.1 3498 14E+00 

9 0.4467487E-19 0.2574701E-01 0.4125847E-19 -0.2574701E-01 

10 0.1002407E-18 0.2574722E+00 0.925939BE-19 -0.2574722E+00 

11 0.10491B5E-18 0.4944374E-01 0.9489564E-19 -0.4944374E-01 

12 0.4149044E-18 0.3449470E+00 0.S497333E-18 -0, 3449470E+00 




C = « 1 

SYSTEM EIGENVALUES 
REAL I HAS INARY 


-0.62898 7 00E-02 
“0.62898700E-02 
“0.3311 771 .56-0.1 
-0.331 17715E-01 
-0.4 7024 58 IE -02 
-0. 47024581 E-02 
-0. 31 19053 7E -01 
-0.31190537E-01 
-0.2 4363 184E-01 
“0. 24363 184E-01 
-0.3 36 2862 IE -03 
-0, 3 362862 9E -03 


0. 1 5666400E+02 
“0. 15666400E+02 
-Q. 6308261 8E+01 
0.63082618E+01 
0. 17995853E+01 
-0. 17995853E+01 
-0, 13384467E+01 
0. 13384467E+01 
-0.94100220E+00 
0. 94100220E+00 
0. 3469791 3E+00 
HU4697913E+00 


DAHPINB 

PERCENT CRITICAL 


0.40148785E-01 
0.40148785E-01 
0. 52498231 E+00 
0, 5249823 lE+00 
0.26130699E+00 
0. 261.3Q699E+00 
0.23297208E+01 
0.23297208E+01 
0,25882004 £+01 
0,2588 20 04 E+01 
0, 96 91 8236 E- 01 
0.96918259E-01 


SYSTEM EIGENVECTORS 


m col l 

REAL IMAGINARY 


1 0. 1000000E+01 0.5055093E-17 

2 -0. 38092 14E “05 -0. 4676142E-07 

3 -0.7158054E+00 -0.4332292E-03 

4 0.9272907E-03 0. 1213413E“04 

5 0.4505361E+00 0.6261089E-02 

6 “0.2257298E+00 -0.3136589E-02 

7 -0. 628987 0E-02 0. 1566640E+02 

8 0.7565426E-06 -0.5967637E-04 

9 0, 1128947E-01 “0. 1121409E+02 

10 -0.1959307E-03 0.1452723E-01 
U -0.1009225E+00 0.7058240E+01 
12 0.5055886E“01 -G.3536354E+01 


SYSTEM EIGENVECTORS 


RON COL 3 

REAL IMAGINARY 


1 “0.7237261 E+00 “0.1663599E-01 

2 -O.354753OE-03 -0.7900585E-05 

3 -0.6962840E+00 -0.2517155E-01 

4 O.1341O51E-01 Q.15O36O8E-03 

5 0. lOOOOOOE+Ol -0.3374579E-16 

6 “0.5065323E+00 -0.71 376 15E-04 

7 -0.8097605E-01 0.4566005E+01 

8 -0.3809035E-04 0.223B137E-O2 

9 -0. 1357294E+00 0.4393176E+01 

10 0.5043897E-03 “0.8460197E-01 

11 -0.3311 77 IE-01 “0.630B262E+01 

12 0. 1632493E-01 0.3195341E+01 


2 

REAL IMAGINARY 


0, lOOOOOOE+Ol -0.2368304E-17 
“0. 38092 14E -05 0.4676142E-07 
-0.7158054E+00 0.4332292E-03 
0.9272907E-03 -0.1213413E-04 
0.4505361 E+00 -0.6261089E-02 
“0.2257298E+00 0.3136589E-02 
-0.6289870E-02 -0.1566640E+02 
0, 7565426E-06 0.5967637E-04 
0.1128947E-01 0.1121409E+02 
-0.1959307E-03 -0.1452723E-01 
“0.1009225E+00 -0.705B240E+01 
0.5055886E-01 0.3536354E+01 


4 

REAL IMAGINARY 


-0,7237261 E+00 0. 1663599E-01 
-0. 3547530E-03 0. 7900585E-05 
-0.6962B40E+00 0.2517155E-01 
0.134105 IE-01 -0.1503608E-03 
0. lOOOOOOE+Ol 0.2648164E-16 
-0.5065323E+00 0.7137615E-04 
“O.8O97605E-O1 -0.4566005E+01 
-0 . 3809035E-04 -0. 22381 37E-02 
-0. 1357294E+00 -0. 43931 76E+01 
0.5043897E-03 0.8460197E-01 
-0.3311 771 E-01 0. 63082 62E+01 
0. 1632493E-01 -0.3195341E+01 
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SYSTEM EIGENVECTORS 

C « .1 (continued) 


RDM COL 5 

REAL IMAGINARY 


1 0.7856006E-01 O. 9365562E-02 

2 -0.8072873E+00 0_ U0$227E-0i 

3 0.2436459E+O9 0.2889396E-01 

4 0, 1000000E+01 0-139591 OE ~ I 6 

5 0. 4 204 28 1 E +00 0.4959143E-01 

6 ~0,43U977EHi0 -0- 2795721E-01 

7 - 0 • 1 7 223 5 5E - 0 1 0,i4133t5£+00 

8 -O, 1605727E-01 **<U452834E+0i 

9 -0.5334288E-0I 0.4383257E+00 

10 -0.47 02 45 BE “02 0.i7?9585E+01 

11 -0,91 221 05E “01 0.7563630E+00 

12 0.5233908E-01 -0,775845*E+0G 


SYSTEM EIGENVECTORS 


RON COL 7 

REAL IMAGINARY 


t 0. 1682317E+O0 -0- i008734E-02 

2 0. 7772059E+00 0. 7777795E-01 

3 0.5502990E+00 -0. 1588694E-02 

4 O, 1563565E+0O 0.3118889E-01 

5 0.1000000E+03 -O. 1791644E-16 

6 -0.7512228E+00 -0.4771140E-01 

7 -0,fc597 37 4E -02 -0.2251377E+00 

8 0, 79860 18E-01 -0.1042675E+01 

9 -0.1929050E-03 -0.7364964E+00 

10 O.10379O2E*QO -0/2U8071E+00 

11 -0.3119054E-01 -0, 1 33844 7E+01 

12 -0. 40428 12E -01 0.1006960E+O1 


SYSTEM EIGENVECTORS 


ROW COL 9 

REAL IMAGINARY 


1 0*1577030E+00 -0.5132384E-03 

2 -0. 668521 7E+00 0.8725714E-01 

3 0.5336765E+00 -0.8123087E-03 

4 -0.7494751E+00 0.6664837E-01 

5 0. 1OOOOO0E+O1 0.1515173E-16 

6 -0. 1702802E+00 -0. 4730163E-01 

7 -0, 43251 06E-02 -0.1483864E+00 

8 0.9839648E-01 0.6269545E+00 

9 -0.1376644E-01 -0.5021710E+00 

10 0.8097586E-01 0.7036339E+00 

11 -0.243631BE-01 -0. 9410022E+00 

12 -0.4036237E-01 0.1613864E+00 


6 

REAL IMAGINARY 


0. 7856006E-01 -0- 9365562E-02 
-O.B072B73E+00 -0. U03227E-01 
0. 2436459E+00 -0 . 2889396E-01 
0. 1000000E+01 0.2748453E-16 
0. 4204281 £+00 -0. 49591 43E-01 
-0,4311 97 7£+00 0.2795721E-01 
-0, 1722355E-01 HU4I3315E+O0 
-0, 1605727E-01 0.1452834E+01 
-0. 531428BE-03 -0. 4383257E+00 
-0.4702458E-02 -0, 1799585E+01 
-0.9122105E-01 -0.7563630E+00 
0.5233908E-01 0.7758456E+00 


8 

REAL IMAGINARY 


0, 1 68231 7E+00 0, 1008734E-02 
0. 7772059E+00 -0.7777795E-O1 
0.5502990E+00 0, 1388694E-02 
0, 1563565E+00 -0.8U8889E-01 
0, lOOOOOOE+Ol 0.220635 1 E- i 6 
-0.7512228E+00 0.4771140E-01 
-0. 6597374E-02 0, 2251377E+00 
0,7986018E-01 0, 1042675E+01 
-0.192905DE-01 0 % 7364964E+00 
0, 1037902E+00 O.2118071E+OO 
-0.3119O54E-O1 0, i 33844 7E+01 
-0.4042812E-03 -O.100696OE+01 


10 

REAL IMAGINARY 


0. 1577030E+00 0.5132384E-03 
-0.6685217E+00 -O.8725714E-03 
O.5336765E+00 0.8123087E-03 
-0. 7494751E+00 -0. 6664837E-01 
O. lOOOOOOE+Ol -0.4171468E-16 
-O.17028O2E+00 0. 4730163E-01 
-0. 4325106E-02 0, 1483864E+00 
0.9839648E-O1 -0.6269545E+00 
-0. 1376644E-01 0.502171QE+00 
0.8097586E-01 -0.7036339E+00 
-0.2436318E-O1 0. 9410022E+00 
-0.4036237E-01 -0. 1613864E+00 
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"SYSTEM EIGENVECTORS 

C = .1 (concluded) 

ROW COL li 12 


REAL IMAGINARY 


t D.2126892E-01 -0. 02472 14E-03 

2 0. 3948017E+00 0. 1367405E-03 

3 0, 7409794E-01 -0. 2873877E-02 

4 O, 742071 4E+00 0. 1848S39E-03 

5 0, 1423528E+00 -0.5522199E-O2 

6 0.10O0OO0E+01 0.1088946E-36 

7 0, 2790087E-03 0. 7380i48E"02 

8 -'0. 1802125E-03 O. 1369879E+00 

9 0,972257 IE -03 O.257 U4OE-01 

10 -O.3O67597E-03 0. 2574832E+00 

11 0.1 86821 6E-Q2 0.4939530E-01 

12 -0.3362862E-03 0.3469791E+OO 


REAL IMAGINARY 


0.2126892E-01 0.8247214E-03 
0. 39480 1 7E+00 -0 . 1367405E-O3 
0.7409794E-01 0. 2873877E-02 
O. 742071 4E+0O -0. 1648839E-03 
0.1423528E+00 0. 55221 99E-02 
0.1OO0O00E+O1 -O. 1133669E-16 
O.2790O87E-O3 -0.7380I48E-02 
-O. 18O2125E-03 -0. 1369879E+00 
0,9722571 E-D3 -0.2571 HOE-Oi 
-0, 3O67597E-03 -0. 2574832E+00 
O. 1868216E-02 -0.4939530E-01 
-0.3362862E-O3 -0. 3469791E+00 
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C = .2 


SYSTEH EIGENVALUES DARPINB 

REAL IMAGINARY PERCENT CRITICAL 


-0.12573562E-01 
-0.12573562E-01 
“Q.66161739E-01 
-0.6616I739E“Ci 
-G.9D785026E-O2 
-0.9O785O26E-O2 
“0.62464502E“01 
-0.62464502E-0! 
“0.49O51767E-01 
-0, 49051 7.67E-01 
-0. 67002982E-03 
HU700299BE-03 


-0, 15666 155E+02 
0.15666155E+02 
0.63060551E+01 
-0.6306«S5tE+0l 
-0.17980485E+01 
0, 17980485E+01 
“O,133143i6E*01 
0, 13314316E+01 
0.94522781E+00 
- 0 . 9452278 1 E+00 
-0.34701547E+00 
0. 34701547E+00 


0.80259375E-01 
0.P0259375E-01 
0. 10491202E+01 
0.1 0491 202E+01 
0.50490222E+00 
0.50490222E+00 
0 » 468637 4 7 E+ 0 1 
0. 468637 47E+01 
0,5iB24388E+0i 
0.51824388E+01 
0, 19308320E+00 
0. 19308324E+00 


BYSTEJI EIGENVECTORS 


RON COL 1 

REAL IMAGINARY 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

it 

12 


SYSTEM EIGENVECTORS 


RON COL 3 

REAL IMAGINARY 


1 -0.7233987E+9O 

2 -0.3549465E-03 

3 -0.6958687E+00 

4 0.134I666E-01 

5 O.iOOOOOOE+Ol 

6 -0.5065353E+00 

7 -0, 1 62023 9E+00 

8 -0.7618740E-04 

9 -Q. 27151 14E+00 

10 0.1008377E-02 

11 -0.66 161 74E-01 

12 0.3261323E-01 


2 

REAL IMAGINARY 


0. 1000000E+01 -0.6860967E-18 
-0.3807563E-05 -0.9348586E-O7 
-0.7157892E+00 -0.8660661E-03 
0.9268317E-03 9.2425712E-04 
0.4502842E+00 0. 1251565E-01 
-0.2256 03 6E +00 -0.6269909E-C2 
-O, 1257356E-01 0. 1566615E+02 
G.I512439E-05 -0.5964870E-04 
0.2256795E-0I -0,1121365E+02 
-0.3916694E-03 0.1451958E-01 
“0. 2017338E+00 0.7054065E+01 
G.1010620E+00 “0.3534263E+01 


4 

REAL IMAGINARY 


“0.72339B7E+00 -0.3328313E-01 
“0. 354946 5E-03 -0, I580565E-04 
-0, 6958687E+00 -0. 5035656E-01 
0. 1341666E-01 0.3006708E-03 
0. lOOOOOOE+Ol 0. 2003064E-16 
“0.5Q65353E+00 -0. 1427249E-03 
-0. 1620239E+00 0.4563994E+01 
-0.7618740E-04 0.2239358E-02 
-0. 27151 14E+00 0.4391 5 18E+01 
O. 1Q08377E-Q2 -0. 84626 ICE-Ol 
-0. 6616174E-01 -0. 6306055E+01 
0.3261323E-01 0.3194249E+01 


0. lOOOOOOE+Ol 
-0.3807563E-05 
-0.7157B92E+00 
0.92683 l 7E-03 
0. 45O2842E+0O 
“0.2256036E+00 
“0. 3257356E-01 
0.1512439E-05 
0.2256795E-01 
“0.3916694E-03 
-0.2017338E+00 
0. 1010620E+00 


0. 92021 66E- 17 
0.93485B6E~07 
0.8660661E-03 
“0, 24257 12£- 04 
-0. 1251565E-01 
0. 62699 09E “02 
-O.1566615E+02 
0«5964870£“04 
0.U21365E+02 
“0.1451958£“01 
“0.7054065E+01 
0.3534263E+01 


0.3328313E-0I 
0.1580565E-04 
O.5935656E-01 
-0.3006708E-03 
0.2848741E-16 
0.1427249E-03 
-0.4563994E+01 
-0.2239358E-02 
-O, 43915 18E+01 
9,846261OE“01 
0.6306055E+01 
“0.3194249E+01 
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SYSTEM EIGENVECTORS 

C = .2 (continued) 


RON COL 5 

REAL IMAGINARY 


1 0.7554 70 7E-01 -0.1810115E-01 

2 -0.8 10530 IE +O0 -0.21462B9E-01 

3 0.2343985E+00 -Q.5586735E-01 

4 0. lOOOOOOE+Ol 0.4626833E-16 

5 0.4O4A4I3E+OO -0. 9592703E-01 

6 -0. 422366 OE +00 O..54110O7E-0I 

7 -0.3323259E-01 -O.135873OE+O0 

8 -0.3123293E-0I 0. 1457567E+01 

9 -O. 1025802E+OO -0.4209528E+00 

10 -0. 9078503E-02 -O.1798O49E+01 

11 -&.1761550E+00 -O.7266939E+O0 

12 0.1011270E+O0 0.7589433E+O0 


SYSTEM EIGENVECTORS 


m COL 7 

REAL IMAGINARY 


1 O. 1679596E+00 -0.2005777E-02 

2 O. 749834 1E+0O O. 1569677E+00 

3 O. 549871 5E+00 -0.3159351E-02 

4 0.1472464E+OO O. 1610129E+00 

5 0. 1 0 0000 OE +0 1 -0.1B41788E-16 

6 -0.7425739E+O0 -0.9525066E-01 

7 -0.131A2O7E-O1 -0. 2235014E+00 

8 8, 1621537E+00 -0.10081 58E+01 

9 -0, 385539 IE-01 -0.7319189E+O0 

10 0.2051 79 9E +00 -Q.2061061E+OO 

11 -0.6246450E-O1 -0, 1331432E+01 

12 -O. 8Q43523E-01 0.994636LE+00 


SYSTEM EIGENVECTORS 


RON COL 9 

REAL IMAGINARY 


1 O. 1577689E+00 O. 1038491E-O2 

2 -0.6469157E+0O -O. 1744151E+00 

3 0. 53378 12E+00 0.1643578E-02 

4 -O. 73357 10E+0O -0. 1334291E+00 

5 0. lOOOOOOE+Ol 0.9812030E-17 

6 -0.1 7895 HE +0O 0.94473O2E-O1 

7 -0.872O455E-02 0. 1490766E+00 

8 O. 1965944E+00 -0.6029273E+00 

9 -0.2773646E-01 0.5044642E+00 

10 0.1621039E+OO -0,686 94 68E+0Q 

11 -0.4905177E-C1 0.9452278E+00 

12 -0.8052066E-O1 -0.1737836E+O0 


6 

REAL IMAGINARY 


0.7554707E-01 O, 1810115E-O1 
-0.81053O1E+0O 0.2146289E-01 
O.2343985E+0O 0.5586735E-OI 
O. lOOOOOOE+Ol 0. 1589034E-17 
0.4O46 4 13E+00 0.9592703E-01 
-O. 4223660E+0O -0.541 1007E-01 
-O, 3323259E -0 1 0.1356730E+00 
-0.3123293E-O1 -O. 1457567E+01 
-O, 1025802E+0O 0.4209528E+O0 
-0.9078503E-O2 0. 1798049E+O1 
-0. 1761550E+O0 O.7266939E+O0 
O.101127OE+00 -O, 7589433E+00 


8 

REAL IMAGINARY 


O.1679596E+O0 0. 2005777E-02 
O. 7499341 E+OO -0.1569677E+00 
0.5498715E+O0 0.3159351E-02 
O.1472464E+0O -0. 1610129E+00 
0, 1 O00OO OE +01 0.2424547E-16 
-O.7425739E+00 0.9525066E-01 
-0. 1316207E-O1 0,223501 4E+ DO 
O, 1621537E+0O 0.1O08158E+O1 
-D, 385539 IE-01 0. 73191S9E+00 
0.2051799E+O0 0.2O6 1061 E+OO 
-0.624645OE-O1 0. 1331432E+01 
-O.8043523E-01 -0.9946361E+OO 


ID 

REAL IMAGINARY 


0. 1577689E+O0 -0. 103849.1E-02 
-O.6469157E+0O 0. 1744151E+00 
0, 53378 12E+00 -0. 1643578E-02 
-O, 73357 lOE+OO 0.1334291E+00 
0. lOOOOOOE+Ol -O. 37188 13E- 16 
-0.17B95UE+00 -O.9447302E-01 
-O.872O455E-02 -0. 1490766E+00 
O.1965944E+0O 0. 6029273E+00 
-O.2773646E-01 -0.5044642E+00 
O.1621O39E+00 O.6868468E+0O 
-0.4905177E-O1 -0, 94522 78E+0O 
-0.8052066E-01 0. 1737836E+00 
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SYSTEM EIGENVECTORS 

C = .2 (concluded) 

ROM COL 11 12 

REAL IMAGINARY REAL IMAGINARY 


1 0. 21 1 7643E-01 0. 1442710E-02 0.2117643E-01 -O. 1442710E-02 

2 0.3948141E+00 -0.2724921E-03 0.3948141E+00 0.2724921E-03 

3 0.7377557E-01 0.57242B0E-02 0.7377557E-O1 -0.5724280E-02 

4 0.7420B89E+00 -0, 3235730E-03 0.7420889E+00 0.3285730E-03 

5 0. 1417332E+00 0.109992BE-01 0.1417332E+00 -0.1099928E-01 

4 0.1000000E+91 0.5258381E-17 0. 1000000E+01 -0. 14401B5E-16 

7 0.555854BE-03 -0.7349449E-02 0.5558568E-03 0.7349449E-02 

8 -0.3590976E-03 -0, 137007 1E+00 -0.3590976E-03 0.1370071E+00 

9 0. 19349B2E-02 -0.2540510E-01 0. 19349B2E-02 0.2560510E-01 

10 -9.6M2416E-03 -0.2575141E+00 -0.4112414E-O3 0.2575141E+00 

11 0.3721953E-02 -0.4919100E-01 0.3721953E-02 0.4919100E-01 

12 -0. 6700299E-03 -0.3470I55E+00 -0.4700299E-03 0.3470155E+00 
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C = .3 


SVSTEH EIGENVALUES 


REAL 


-0, 18845095E-01 
-0. 18845O95E-01 
-0. 990583 i8E~01 
-0.990583I8E-01 
-0.12855629E-01 
-0.12855A29E-01 
-0.93807180E-01 
-Q. 93807 180E-01 
-0.74435443E-01 
-9.74435443E-01 
-0.99874304E-03 
-0.99874328E-03 


IMAGINARY 


0.1 5665746E+02 
“0, 15665746E+92 
0. 6302381 8E+01 
-0. 630238 18E+0 l 
0.17956874E+01 
-0.17956874E+01 
-0.13189978E+01 
0.13189978E+01 
0.95275050E+00 
-O.95275O50E+90 
0.347075 48E+ 00 
-0.34707548E+00 


DAMPING 

PERCENT CRITICAL 


0, 12029483E+00 
0. 12029483E+00 
0. 15715659E+01 
0.15715659E+01 
0.7158 9855E+ 00 
0.71589855E+00 
0.70940857E+01 
0. 70940857E+01 
Q.77B89549E+01 
O,77889549E+0i 
0.28775842E+00 
0.28775849E+G0 


SYSTEM EIGENVECTORS 


ROM COL 1 

REAL IMAGINARY 


1 O.IOOOOOOE+Ol -0.9581637E-17 

2 -0.38O4816E-O5 -0. 1401363E-06 

3 -0.7157621E+00 -0. 1298126E-02 

4 0.9260674E-03 0.3635784E-04 

5 M49B650E+OO 0, 1875719E-01 

6 -0.2253936E+00 -0, 9394705E-02 

7 -0.1B84510E-01 O. 1566575E+02 

8 O.2267O42E-05 -0.5960264E-04 

9 0,338247 IE -01 -O. 1121292E+02 

10 -0.5870246E-03 0.1450685E-01 

11 -0, 302323 1E+00 0. 70471 17E+01 

12 0. 1514540E+00 -0.3530782E+01 


SYSTEM EIGENVECTORS 


ROM COL 3 

REAL IMAGINARY 


1 -0.7228538E+00 0.4995251E-01 

2 -0.3552705E-03 0.2371960E-04 

3 -0.6951773E+00 0.7556851E-01 

4 0.1342694E-01 -O.450878OE-O3 

5 0. iOOOOOOE+Ol 0.1 22 1 083E- 1 6 

6 -O.5G65403E+OO 0.2140181E-03 

7 -0.2432151E+00 -0.4560649E+01 

8 -0. 1 142975E-03 -0.2241400E-02 

9 -0.4073985E+00 -0. 4388759E+01 

10 0.1511556E-02 0.8466636E-01 

11 -0.9905832E-01 0.6302382E+0! 

12 0.4882821E-01 -Q.3192432E+01 


REAL IMAGINARY 


0, IOOOOOOE+Ol 0.30981 08E- 16 
-Ot 380481 6E-05 0.1401363E-06 
-0.71 5762 IE +00 0. 1298126E-02 
0.9260674E-03 -0.3635784E-04 
0.4498650E+00 -0.1875719E-01 
-0.2253936E+00 0.9396705E-02 
-0,1 8845 10E-01 -0. I566575E+02 
0.2267042E-05 0.5960264E-04 
0, 3382 471 E -01 0.1121292E+02 
-0.5870246E-03 -0.1450685E-01 
-0. 3023231 E+00 -0.7047117E+01 
0.1514540E+00 0.3530782E+01 


4 

REAL IMAGINARY 


-0. 7228538E+00 -0.4995251E-01 
-0.3552705E-03 -0.237I960E-04 
-0.6951773E+00 -0.7556851E-01 
0.1342694E-01 0.4508780E-03 
0, IOOOOOOE+Ol -0.3480289E-16 
-0.5065403E+00 -0.2140181E-03 
-0.24321 5 1E+00 0.4560649E+01 
-0.1142975E-G3 0.2241400E-02 
-0.4073985E+00 0.4388759E+01 
0.1511556E-02 -0.8466636E-01 
-0. 9905832E-01 -0.6302382E+01 
0.4882821E-01 0. 3192432E+01 
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■SYSTEM EIGENVECTORS 


C = .3 (continued) 


RON COL 5 

REAL IMAGINARY 


1 0.7090664E-01 0.2567787E-Q1 

2 -0.8 156 1 48E+00 0.3075682E-01 

3 O.22O14Q6E+O0 0.7930222E-01 

4 0, 1000000E+01 -0.2225325E-16 

5 0 . 38027288 +00 (U362540E+D0 

6 -0 ,40871 31 E +00 -0. 7692A20E-01 

7 -0. 4702098E-01 0. 1269960E+00 

8 -0. 4474439E-01 -0. 1464985E+01 

9 -0, 1452320E+00 0.3942842E+00 
10 -0. 1285563E-01 0.i795687E+01 
1! -O.2495582E+0O 0.6810995E+00 
12 0. 1433897E+00 -0. 73293 19E+00 


SYSTEM EIGENVECTORS 


m COL 7 

REAL IMAGINARY 


1 0. 1674B67E+00 -0.2974251E-02 

2 0. 7033490E+00 0. 23S8106E+00 

3 0, 5491282E+00 -0.468580 IE-02 

4 O. 1301355E+00 0.2383O35E+OO 

5 0. iOOOOOOE+Ol 0.3084555E-16 

6 -0.7273081E+00 -0.1424926E+00 

7 -0.1963448E-01 -0.2206355E+00 

8 0,2490 11 5E +00 -0. 95011 80E+00 

9 -0.57 6927 3E-01 -0.7238593E+00 

10 0.3021142E+00 -O.194O03OE+O0 

11 -0. 938071 8E -01 -0, 13I8998E+01 

12 -0.1197207E+00 0.9726846E+00 


SYSTEM EIGENVECTORS 


RON COL 9 

REAL I HAS INARY 


1 0. 1578884E+00 0.1589878E-02 

2 -O. 609268 1E+0O -0.2613393E+00 

3 0.5339708E+00 0.2516092E-02 

4 -0.7059264E+00 -0.2004835E+00 

5 0. 1000000E+01 -0.1B07907E-16 

6 -0.1941389E+00 0. 141.3734E+00 

7 -0.1326725E-01 0.15D3099E+00 

8 0.2943423E+00 -0.5610276E+00 

9 -0.4214356E-01 0.5085537E+00 

10 O.2435567E+O0 -0.6576487E+00 
U -0.7443544E-01 0.9527505E+00 
12 -0.1202428E+00 -0.1 95489 1E+00 


6 

REAL IMAGINARY 


0.7090664E-01 -0.25677S7E-01 
-0. 8156148E+00 -0. 3075682E-01 
0.220H06E+OD -0.7930222E-01 
0. lOOOOOOE+01 -0.19298B0E-16 
0.3802728E+00 -O. 1362540E+00 
-0,40871 3 lE+00 0.7692620E-01 
-0.4702098E-01 -0.1269960E+00 
-0.4474439E-01 0.1464985E+01 
-O.1452320E+80 -0.3942B42E+00 
-0.1285563E-01 -0. 1795687E+QI 
-0, 2495582E+00 -0.681C995E+00 
0.1433897E+00 0.7329319E+00 


8 

REAL IMAGINARY 


0, 1674867E+00 0.2974251E-02 
0. 7033490E+00 -Q. 23881 06E+00 
0.5491282E+00 0.4685801E-02 
0. 1301355E+00 -0.2383035E+OO 
0. lOOOOOOE+Ol -0.6U2190E-16 
-0. 7273081E+00 O. 1424926E+00 
-0.1963448E-01 0.2206355E+00 
0, 24901 15E+00 0.9501 180E+00 
-0.5769273E-01 0.723B593E+00 
O. 30211 42E+00 0. 1940030E+00 
-0, 93807 18E -01 0. 131899BE+01 
-0.1197207E+00 -0.9726846E+00 


10 

REAL IMAGINARY 


0, 1578884E+00 -0. 1589878E-02 
-0.60 9268 IE +00 0.2613393E+00 
0. 533970BE+00 -0. 251 6092E-02 
-0.7059264E+00 0.2004835E+00 
0.1000000E+01 -0.2429290E-17 
-0.1941389E+00 -0. 1413734E+00 
-0. 1326725E-01 -0.1503099E+00 
0.2943423E+00 0.5610276E+00 
-0.4214356E-01 -0.5085537E+00 
0.2435567E+00 0.6576487E+00 
-0.7443544E-01 -0.9527505E+00 
-0. 1202428E+00 0. 1954891E+00 
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SYSTEM EIGENVECTORS 

C = .3 (concluded) 

ROW COL 11 12 

REAL IMAGINARY REAL IMAGINARY 


1 0.21G2385E-01 

2 0.39484Q0E+00 

3 0,73243766-01 

4 0.7421177E+00 

5 0, 14071 12E+00 

6 0, 1000000E+01 

7 0.828432OE-O3 

8 -0,53535576-03 

9 0.2886809E-02 

10 “0,91 121 62E -03 

11 0,55470536-02 

12 -0.9987432E-03 


-0, 244739 IE-02 
0.4062861E-03 
-0.8S28292E-02 
0, 48989 68E-03 
-0, 1 63871 8E-01 
-0.IO43545E-16 
0.7299306E-02 
0. 1370389E+00 
0.2542-9636-01 
0.25757 04 E+OQ 
0.4 88537 6E- 01 
0.3470755E+00 


0.2102385E-01 
0. 3948400E+00 
0.7324376E-01 
0, 7421 177E+00 
0.14071126+00 
0. 10QOQ0OE+O1 
0. 828432 OE -03 
-0.5353557E-O3 
0. 2886 80 9E “02 
-0.9112162E-03 
0.5547053E-02 
-0.9987432E-03 


0.2447391E-02 
-0.4062861E-03 
0.B528292E-02 
-0.4898968E-03 
0.1638718E-01 
“0.1699487E-16 
“0.7299306E-02 
“0. 13703B9E+00 
-O.2542963E-01 
“0.2575704E+00 
“0. 488537 6E“ 01 
-0.3470755E+00 
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C = .4 


SYSTEM EIGENVALUES DAMPING 

REAL I HAS I NARY PERCENT CRITICAL 


-0.25097802E-01 
-0.25097802E-01 
H0.13J73314E+00 
-0* 131 73.31 4E tOO 
-O.I5858107E-01 
-0.15858I07E-01 
-0.124W314E+00 
-0.12499314E+O0 
~O.10O99797E+00 
-0.IO099797E+00 
~0. 132OO455E-02 
-0/132O0459E-02 


-O.15.W51.74E+02 
O. 15665 174E+02 
0.6297249OE+01 
-0. 6297249OE+01 
0. 17927596E+01 
-0. 17927596E+01 
-0, 12998665E+01 
O.12998665E+01 
O.96449645E+D0 
“0.96449645E+00 
*0,34715832E*0O 
O.34715832E+O0 


0, 1602138OE+0O 
O. 16021380E+QO 
O.20914582E+O1 
O.20914582E+O1 
0.88452947E+0O 
0.88452946E+O0 
0.95716937E+01 
0.95716937E+01 
0.10414631E+02 
O, 10414631E+02 
0.38024O35E+0O 
O.38O24044E+00 


SYSTEM EIGENVECTORS 


RON COL 1 2 

REAL IMAGINARY REAL IMAGINARY 


1 0.100DO0OE+01 O.9025983E-17 

2 ■"0.3800975E-05 O. I866764E-06 

3 -0.7157242E*0O O.1729O0OE-02 

4 0.9249993E-03 -0.4842527E-O4 

5 0.4492792E+O0 -0. 2497922E-01 

6 "O.2251002E+00 O.1251373E-01 

7 -0. 250978OE-01 -O. 1566517E+02 

8 0.3019714E-05 0.5953825E-O4 

9 O.4504819E-01 0. 1 121 190E+02 

10 -0.7818O57E-03 -0.I4489O6E-01 

11 -0, 402579 7E +0O -0.7037410E+01 

12 0.2016792E+0O 0.3525920E+OI 


SYSTEM EIGENVECTORS 


RON COL 3 

REAL IMAGINARY 


1 -0. 7220926E+OQ 0.6665542E-01 

2 -O.3557273E-03 0.3164666E-O4 

3 -0.69421 1OE+00 0.1008209E+00 

4 0.I344137E-01 *0.60092511*03 

5 0.10O0O00E+01 -0.992O450E-17 

6 -Q.5065473E+OO 0.2852252E-03 

7 -0.3246223E+00 -0.4555977E+01 

8 -0, 1524258E-03 -O.2244272E-02 

9 -O. 5434439E+00 -0.4384901E+01 

10 O.2013502E-02 0.8472279E-01 

11 -O.1317331E+O0 0.6297249E+O1 

12 0.6493293E-01 -0.3189892E+O1 


0.10O0000E+01 0. 1409463E-16 
-0.3800975E-05 -0.1866764E-06 
-0.7157242E+00 -0. 1729000E-02 
0.9249993E-03 O.4842527E-04 
O.4492792E-V0O 0.2497922E-01 
-0.2251OO2E+00 -0.1251373E-01 
-0.2509780E-01 0, 1566517E+02 
0. 3019714E-O5 -O. 5953825E-04 
0. 450481 9E -01 *0,1121 19OE+02 
-0.7818057E-O3 0.1448906E-O1 
-0. 4025797E+00 0.7037410E+01 
0.2016792E+00 -0.352592OE+01 


4 

REAL IMAGINARY 


*0, 7220926E+O0 -0.6665542E-O1 
*0. 3557273E-03 -0.3164666E-04 
*0. 69421 10E+0O -O.1008209E+00 
0, 13441 37E-01 0.6009251E*03 
0.1000000E+01 0.55511 15E-16 
*0,5065473E+00 ~0.2852252E*03 
-0.3246223E+00 0.4555977E+01 
-0.1524258E-O3 0.2244272E-02 
-0.5434439E+O0 0.4384901E+01 
0.2013502E-02 -0.8472279E-O1 
-O.1317331E+00 -O.6297249E+01 
O.6493293E-01 0.3189892E+O1 
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SYSTEM EIGENVECTORS 

C = .4 (continued) 


RON COL 5 

REAL IMAGINARY 


1 0 , 65 1 3201 E -01 0.3174912E-OI 

2 -0.822O986E+OO 0.3851264E-01 

3 0,2023 7 13E+00 0,98i28?9£-0i 

4 0, 1 OOOOOOE+O 1 -0.6071532E-17 

5 0, 3498563E+O0 0, 1687363E+00 

6 “0.3916367E+OO -0.9537219E-01 

7 -0 ,579514 OE -0 i 0, 1162626E+D0 

8 “0.5600698E-01 -0, 1474436E+01 

9 ”0.I791309E+O0 0, 3612469E+00 

10 -0. 1 585811 E-01 0. 1792760E+01 

11 -0.30805 l 6E+00 C.6245324E+00 

12 0, 17719OOE+0O -O.7Q05981E+OO 


SYSTEM EIGENVECTORS 


RON COL 7 

REAL IMAGINARY 


1 0. 1667B1 6E+00 -0.3886295E-02 

2 0 . 6 3 56043E + 00 0.3242755E+00 

3 0.5480 3 95E+00 -0. 6124582E-02 

4 0, 1018152E+OO 0.3122430E+00 

5 0.1 OOOOOOE+O 1 Q, 59902 17E” 16 

6 -0. 7038785E+00 -0, 1894088E+00 

7 -0.2589823E-01 -0.2163081E+00 

8 0. 3420687E+00 -O.8667329E+O0 

9 -0. 76459B2E-0I -0.7U5867E+00 

10 0.3931 481E+00 -0.1713744E+00 
31 '0.1 249931 £+00 -0.1299867E+01 
12 ”0. 158226 2E+00 0.9386229E+00 


SYSTEM EIGENVECTORS 


RON COL 9 

REAL IMAGINARY 


1 0, 1580804E+00 Q.21870Q2E-02 

2 ”0.55243 L5E +00 -0.3477948E+00 

3 0.5342756E+00 0.3460739E-02 

4 -0.6643558E+00 -0.2679729E+00 

5 O.IODOOOOE+Ql -0. 1409463E”17 

6 -0.2172453E+00 0.1878309E+00 

7 ”0.1807516E”01 0. 1522471E+00 

8 0.391 24 13E+00 -0.4976917E+00 

9 -0.5729862E-01 0.5149574E+00 

10 0.3255575E+00 -0.6137041E+00 

11 -0.1009980E+00 0.9644965E+00 

12 ”0. 1592209E+00 -0.2285028E+00 


6 

REAL IMAGINARY 


0.6513201E-01 -0.3174912E-01 
“0.8220986E+00 -0.3851264E-01 
0.2023713E+00 -0.9812899E-01 
0. I000000E+01 ”0.1614445£“17 
0. 3498563E+00 -0, 1687363E+00 
-0.3916367E+00 0.9537219E-O1 
-0.5795140E-01 -0. U62626E+00 
-0.5600698E-01 0.1474436E+01 
“0. 1791309E+0D -0.3612469E+00 
-0, 1585811E-01 “0. 1792760E+01 
-0.3080516E+00 -0.6245324E+00 
0, I771900E+00 0. 700598 1E+00 


8 

REAL IMAGINARY 


0.1667816E+00 0.3886295E-02 
0.6356043E+Q0 -0.3242755E+00 
0. 54801 95E+00 0.6124582E-02 
0.101 81 52E+00 “0.3122430E+00 
0. 1000000E+01 ”0. 1687290E-17 
“0.7038785E+00 0. 1894088E+00 
“0, 2589823E-01 0.2163081E+00 
0, 3420687E+00 0.8667329E+00 
-O, 7645982E-01 0.7115867E+00 
0.3931 481 E+00 0. 1713744E+00 
”0,1249931 E+00 0.1299867E+01 
”0. 1582262E+00 -Q.9386229E+0Q 


10 

REAL IMAGINARY 


0.1580804E+QO -0.2187002E-02 
-0.5524315E+00 0.3477948E+00 
0.5342756E+00 -0.3460739E-02 
“0.6643558E+00 0.2679729E+00 
0, 1000000E+01 0. 42094 15E” 16 
-0.2172453E+00 -0.1878309E+00 
-0.1B07516E-01 -0 . 152247 1E+00 
0,391 241 3E+00 0.4976917E+00 
-0.5729862E-01 -0.5149574E+00 
0.3255575E+00 0.6137041E+00 
”0.1009980£+00 -0.9644965E+00 
”0. 1592209E+00 0.2285028E+00 
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SYSTEM EIGENVECTORS 

C = ,4 (concluded) 


ROW COL 11 


REAL IMAGINARY REAL 


1 0, 2081346E-01 0.3232503E-02 

2 0.3948730E+00 -0.5371934E-03 

3 0.7251047E-01 0.1126408E-01 

4 0.7421576E+00 -0.6477329E-03 

5 0.1393019E+00 0.2U4394E-01 

6 O.IOOOOOOE+Ol -0.46918B5E-16 

7 0.1094716E-G2 -0.7.229833E-02 


0.2081346E-01 
0.3948730E+00 
0.7251047E-01 
0.7421576E+00 
0.1 3930 19E+00 
O.IOOOOOOE+Ol 
0, 1094716E-02 


8 -Q.70774UE-03 -0.1370827E+OQ -0.7077416E-03 

9 0.3814703E-02 -0.251874BE-01 0.3814703E-02 


10 ”0. 1204548E-02 -Q.2576453E+00 -0, 1204.548E-02 


M 0 , 7329990E-02 -0 . 4838838E-0 1 0.7329990E-02 
12 -0,1 32004 6E-02 -0, 3471 583E+00 -0. 132G046E-02 


12 

IMAGINARY 


-0.3232503E-02 
0, 5371934E-03 
”0, 1126408E-01 
Q.6477329E-03 
-0.21 64394 E-0 1 
0, 27891 10E- 16 
0.7229833E-02 
0.1 370827 E ^-00 
0.25I8748E-0.1 
0.2576453E+00 
0. 483B838E-01 
0.3471583E+00 
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C = .5 


SYSTEM EIGENVALUES 

DAMPING 

REAL IMAGINARY 

PERCENT CRITICAL 


~0. 31325969E-01 

-0. 15664441E+02 

0. 19998 lOlE+OO 

-0.31325969E-01 

0. 1 566444 1E+02 

0.1 9998101 E+00 

-0. 16411 189E+00 

-0.62906666E+01 

0.26079283E+01 

-0.1641 13 89E+00 

0. 62906666E+01 

0.26079283E+01 

-0. 1 802697 OE -01 

0.17895488E+01 

0,1 007296 IE+01 

-O.1802697OE-01 

-0.17895488E+01 

0.1007296 IE+01 

-0. 15540295E+00 

-0. 12716557E+01 

0. 12130279E+02 

-0.15540295E+00 

0.12716557E+01 

0.12130279E+02 

-0. 12950050E+00 

-0. 98242523E+00 

0, 33068666E+02 

-0.12950050E+00 

O.98242523E+0O 

0. 33068666E+02 

-0, 1 6317:1 35E-02 

0.34726287E+00 

0.46987331 E+00 

-0.16317139E-02 

-0, 34 726287 E+00 

0.46987343E+DO 


SYSTEM EIGENVECTORS 


ROW COL 1 

REAL IMAGINARY 


1 O.iOOOOQOE+Ol 0.9730714E-17 

2 -0.3796046E-05 0.233G692E-O6 

3 -0. 7156756E+00 0. 215831 5E-02 

4 0.9236290E-03 -0.6O44840E-O4 

5 0.4485280E+00 -0.3117533E-01 

6 -0.2247239E+00 0.1563777E-01 

7 -0, 31 32597 E-01 -0. 1566444E+02 

8 0. 37698 14E -05 0.5945563E-04 

9 0.5622803E-01 0- 1121059E+02 

10 -0.9758239E-03 -0.1446624E-01 

11 -0, 50239 4 7E +00 -0, 7024963E+01 

12 O.2516834E+0O 0.3519685E+01 


SYSTEM EIGENVECTORS 


ROW COL 3 

REAL IMAGINARY 


1 -0.721 11 66E +00 -O.B340300E-01 

2 -0.3563202E-03 -0.3959088E-O4 

3 -0 , 69297 15E +00 -0.1261274E+00 

4 0. 1345998E-01 0.7507509E-03 

5 0. 1000000E+01 0.2379824E-16 

6 -O.5O65563E+0O -0. 35631 38E-03 

7 -0.4063166E+00 0.4549991E+01 

8 -0.3905767E-03 0.2247989E-02 

9 -0. 6797008E+00 0.437?952E+0i 

10 0 . 25 1 3780E-02 -0.8479547E-01 

11 -0, 16411 19E+00 -O, 629066 7E+01 

12 0.8089046E-01 0.3186635E+01 


2 

REAL IMAGINARY 


0. 1000000E+01 0.8565I97E-17 
-0 . 3796046E-05 -0.2330692E-06 
-0.7156756E+00 -0.21 583 15E-02 
Q.7236290E-03 0.6044840E-04 
0.4485280E+00 0.3117533E-01 
-0.2247239E+00 HU561777E-01 
-0.3132597E-01 0. 1566444E+02 
0.37698I4E-05 -0.5945563E-04 
0.5622803E-01 -0.U21059E+02 
-0.9758239E-03 0. 14466 24E-01 
-0. 5023947E+00 0.7024963E+01 
0.2516834E+00 -0.3519685E+01 


4 

REAL IMAGINARY 


-0,721 11 66E+00 0.S340300E-01 
-0. 3563202E-03 0.3959088E-04 
-0.6929715E+00 0.1261274E+00 
0. 1345998E-01 -0.7507509E-03 
0, lOOuOOOE+Ol 0.2219904E-16 
-0. 5O65563E+00 0.3563138E-03 
-O.4O63166E+O0 -0.454 999 iE+01 
-0.1905767E-03 -0.2247989E-02 
-0 , 6797008E+00 -0. 43799526*0 1 
0.2513780E-02 0.8479547E-01 
-0. 16411 I9E+O0 0. 629066 7E+ 01 
0. 8089046E-01 -0. 3186635E+0I 
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SYSTEM EIGENVECTORS 

C = .5 (continued) 

m COL 5 

REAL IMAGINARY REAL 


1 O,5877046E-Oi 0, 34I8999E-01 O,5377046E-O! 

2 -0,82944685+00 0.4451666E-OL -0,829446BE+0O 

3 0. 1827622E+0Q 0.11195035+00 0,18276225+00 

4 0. lOOOOOOE+Ol -0,21616235-17 O.lOQOOOOE+Ol 

5 Q,3162310E+00 O.1926707E+OO 0,31623i0E+Q0 

6 ~0,3727132E+00 -0. 1090369E+00 -0.3727132E+00 

7 -0.658232 IE -01 0. L045202E+00 -0,65823215-01 

8 -0.6 4 71 23 2E - 0 1 -0, 1485138E+D1 -0.6471232E-01 

9 -0.20363525+00 0,325Q437E+00 -O.2O363525+O0 

10 -0.1802697E-01 0, 17895495+01 -0.18026975-01 

11 -0.3504943E+00 0.5624378E+0O -0.3504943E+00 

12 0 . 20 18 45 BE ■*- 0 C -0.66502285+00 0.2018458E+00 


SYSTEM EIGENVECTORS 


RON COL 7 


REAL IMAGINARY REAL 


1 0. 1657902E+OO -0, 4693972E-02 0.16579025+00 

2 0.5410750E+00 O.4141615E+O0 0.5410758E+00 

3 0. 5464592E+00' -0. 74007045^02 0.5464592E+00 

4 0.5654809E-01 0.3824349E+00 O.5654809E-O1 

5 O. lOOOOQOE+Ol -0,14772255-17 0,10000005+01 

6 -0.66921125+00 -0.23614465+00 -0.66921125+00 

7 -0.31733405-01 -0. 21009865+00 -0. 31733405-01 

8 0.4425861E+00 -0.7524240E+00 0.4425861E+00 

9 -0.9433252E-0I -0.6937578E+00 -0.94332525-01 

10 0,47753775+00 -0.13134125+00 0.47753775+00 

11 -0.15540305+00 -0,12716565+01 -0.1554030E+00 

12 -0.19629725+00 0,88770375+00 -0. 19629725+00 


SYSTEM EIGENVECTORS 


RON COL 9 


REAL IMG I NARY REAL 


1 0,15838635+00 -0.2862923E-02 0,15838635+00 

2 -0.4697707E+00 0.4332587E+00 -0,46977075+00 

3 0,53476095+00 -0,4529626E-02 0.5347609E+00 

4 -0.60425845+00 0.3360869E+00 -0,60425845+00 

5 0, i OOOOOOE+O i 0,35236575-17 0,10000005+01 

6 -0.251 19UE+00 -0.233579DE+00 -0.25119UE+00 

7 -0. 23323725-01 -0,15523205+00 -0.23323725-01 

8 0,4864798E+00 M054074E+OO 0.4864798E+00 

9 -0,73701825-01 -0.524 77605+00 -0,73701825-01 

10 0, 4084320E+00 0,55011535+00 0,40843205+00 

11 -0. 1295005E+00 -0.9824252E+00 -0.1295005E+00 

12 -0. 1969445E+00 0.27702515+00, -.0,19694455+00 


IMAGINARY 


•0,36189995-01 
•0,44516665-01 
■0,11195035+00 
0.3686287E-17 
■0. 1926707E+00 
0. 10903695+00 
•0,10452025+00 
0,14851385+01 
•0.32504375+00 
■0.1789549E+01 
•0.56243765+00 
0.66502285+00 


8 

I HA BI NARY 


0,4693972E-02 
-0,41416155+00 
0.7400704E-02 
-0.3B24349E+00 
0, 252348 IE-1 6 
0.23614465+00 
0,21009865+00 
0,75242405+00 
0 .69375785+00 
0, 13134 12E+00 
0.12716565^1 
-0.8877037E+00 


10 

IMAGINARY 


0, 2862923E-D2 
-0.4332587E+00 
0, 4529626E-02 
-0.3360869E+O0 
0.6209768E-16 
0,23357905+00 
0. 1552320E+00 
-0,40540745+00 
0,52477605+00 
-0.55011535+00 
0.98242525+00 
-0, 277025 1E+00 
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SYSTEM EIGENVECTORS 

C = .5 (continued) 


m col ii 

REAL IMAGINARY 


12 


REAL 


IMAGINARY 


1 O,2O54838E-0i -0. 3992236E-02 

2 0,39491 47E+0O 0,6643424E-03 

3 0, 7i 58656E-01 -0. 1391 I40E-01 

4 0, 742207 9E+00 U.S010287E-03 

5 '0. 1375263E+00 -0.2673064E~Oi 

6 0. 1-090000E+01 -0.37865 76 E~ 1 6 

7 0. 1352826E-02 0.7142202E-02 
•8 -O..B75O893E-03 0. 1371381E+00 

9 0.47141 0 4£ -02 0.24fl8205£-0i 

10 -0, 14S9233E-02 0.25773WE+00 

11 0. 9058 I56E-02 0. 47001 38E-01 
.12 “0, 16317 14E-02 0.3472629E+00 


0.2054838E-01 0.3992236E 02 
0 , 394 9147E +00 -0. 6643424 E-03 
0, 7158656E-01 0,1391 140E-01 
DJ422079E+00 -0.801 0287 E- 03 
0. 1375263En)0 0, 2673064E-01 
0, lOOOOOOE+Ol O, 36429 19E- 16 
0, 1352826E-02 -0.7142202E-02 
- 0 , 8750 8 9 3E - 03 HU371381E*00 
0. 47141 04E-02 -0.2488205E-0! 
-0, 1489238E-02 -0, 2577399E+00 
0.9058 156E-02 -0.4700138E-01 
-0, 1631714E-02 -0. 3472629E+00 
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C = 1.0 


system eigenvalues damping 

REAL IMAGINARY PERCENT CRITICAL 


-0, .6 1.89.2 4 V7E-0 1 
-0.61892497E-01 
-Q. 3S892409E+00 
-O.31892409E+0O 
-0. 19 751 5S0£ “01 
-0.19751551E-01 
-0, 51 26347 3E+00 
-0*5 326347 3£ +00 
-0.83S13675E-01 
-0.83813675E-01 
-0. 29834 595E-02 
-0.2983460 IE -02 


-0, 15658408E+02 
0, 156584 08£+ 02 
-0, 62365954E+01 
O.62365954E+01 
0.17758260E+01 
-0, 17758260E+01 
-O.10174O53E+O1 
0, 101 7 405 3E +01 
0, 1 12626 ISE+Oi 
-0.1 12626 15E+01 
-O.348O5529E+0O 
0.34805529E+00 


0.39526378E+00 
O.39526378E+0O 
0.51070798E+01 
0. 51070798E+01 
0. 11121 770E+01 
0.H121771E+01 
0, 4499726 1E+02 
0. 44997261 E+02 
0. 74212386E+01 
0.74212386E+01 
0, 857 1482 6 E+ 00 
0,8571 48 46E+00 


SYSTEM El SEN VECTORS 


row COL 1 


2 


REAL 


IMAGINARY REAL IMAGINARY 


1 O. lOOOOOOE+Ol 0.207523 IE- 17 

2 -0.3755374E-05 0.4615726E-06 

3 -0.7152753E+O0 0.4268330E-02 

4 0.9123443E-03 -O. 119525 IE-03 

5 0,4 4235 18E+00 -0.6154943E-01 

6 -0.2216302E+00 0.3083425E-01 

7 -0.6189250E-01 -O. 156584 1E+02 

8 0, 7459922E-05 0.5877461E-04 

9 O.lll 1054E+00 O. 11 19981 E+02 

10 -0. 1928039E-02 -0. 1427846E-01 

11 -0.9911443E+OO -0.6922715E+Q1 

12 0.4965324E+00 0.3468468E+01 


0, lOOOOOOE+Ol -0. 14924 72E- 17 
-0. 375537 4E -05 -0.4615726E-06 
-0.7152753E+00 -0.426B330E-02 
0.9123443E-03 0.1195251E-03 
0, 442351 8E +00 0.6154943E-01 
-0, 22163O2E+O0 -0.3083425E-01 
-0.6189250E-01 0. 1565841E+02 
O. 7459922E-05 -0.587746 IE -04 
0. 1111054E+0O -O.l 119981 E+02 
-0. 1928039E-02 0. 1427846E-01 
-0,9911 44 3E+00 0. 69227 15E+01 
O. 4965324E+00 -0, 3468468E+01 


SYSTEM EIGENVECTORS 


RON' COL 


4 


REAL IMAGINARY REAL 


IMAGINARY 


1 -0.7131O97E+00 

2 -O.3635O5OE-03 

3 -0.6827750E+00 

4 0.33&1798E-01 

5 0. lOOOOOOE+Ol 

6 -0.5066327E+00 

7 -0.8216029E+00 

8 -0.3816612E-03 

9 -0. 1366072E+01 

10 0.49 7288 IE-02 

11 -0. 318924 1E+00 

12 0. 1571 586E+00 


-O. 1 68205 7E+0O 
-0.7968353E-04 
-0.2539568E+00 
0.1493761E-02 
0.3699840E-16 
-O.708528OE-O3 
0.45O1O22E+O1 
0.227 997 4E -02 
0. 43391 85E+01 
-0.8540625E-OI 
-0.6236595E+01 
0.3159889E+01 


-0.7131097E+00 
-0.3615050E-O3 
-0.6827750E+00 
0.1361798E-01 
0, lOOOOOOE+Ol 
-0.5066327E+00 
-0.8216029E+00 
-0.38 166 12E-03 
-0. 1366072E+01 
0.49 7288 IE -02 
-0.31 8924 IE +00 
O. 1571 586E+00 


0. 1682057E+00 
0.7968353E-04 
0.2539568E+00 
-0.1493761 £-02 
0. 1 79571 OEM 7 
0.7085280E-03 
-0.4501O22E+O1 
-0.2279974E-02 
-O,4339i85E+0i 
0, 854O625E-01 
0.6236595E+01 
-0.3159889E+01 
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C - 1.0 (continued) 


6 


ROW COL 5 

REAL IMAGINARY REAL 


1 O,31i5362E-0i 0.4020393E-01 0.3115362E-01 

2 -0.8639792E+00 0.5255843E-01 -0.8639792E+00 

3 0.9723323E-01 0. 1248190E+00 0.9723323E-O1 

4 0, 1000000E+01 -0.271 0505E- 18 0. lOOOOOOE+Oi 

5 O.168B589E+O0 0.2156128E+00 0.1688589E+00 

6 -0.2891885E+00 -0. 1227091E+00 -0.2891885E+00 

7 -0.7201052E-01 0.5452933E-Q1 -0.7201052E-01 

8 -0. 7626970E-01 -0.1535315E+01 -0.7626970E-01 

9 ~0, 223577 4E+00 0, 1702039E+00 -O.2235774EH10 

10 HU975155E-01 0= 1775826E+01 -0. 1975155E-01 

11 -0, 3862261 E+00 9.2956053E+00 -0.3862261E+00 

12 0.2236219E+00 -O.5111248E+O0 0.2236219E+00 


SYSTEM EIGENVECTORS 


RON COL 7 

REAL I HAGINARY REAL 


1 0, 1559523E+00 -0. 1154731E-QI 0.1559523E+0Q 

2 -0.1951235E-01 0-1781620E+00 -0. 1951235E-01 

3 0.5309572E+00 -0. 1828979E-01 0.5309572E+00 

4 -0.20979 83E+00 0. 1983724E+00 -0. 2097983E+00 

5 0. 1000000E+01 0.1219727E-16 0. 1000000E+01 

6 -0.4449782E+00 -0.1534676E+00 -0. 4449782E+00 

7 “0.9169487E-01 -0. 1527472E+00 -0, 9169487E-01 

8 0. 1912657E+00 -0.714BG06E-01 0. 1912657E+00 

9 -0.2907952E+00 -0.5308227E+00 -0.2907952E*00 

10 0. 309375 1E+00 0.U17573E+G0 0.3G93751E+00 

11 “0.5126347E+00 -0.1017405E+01 -0.5i26347E+00 

12 0.7197249E-01 0.5313960E+00 0.7197249E-01 


SYSTEM EIGENVECTORS 


RON COL 9 

REAL IMAGINARY REAL 


1 D, 14047 19E-02 0.1021558E+00 0. 1404719E-02 

2 0.1000000E+01 0.2385245E-17 0. 1000000E+01 

3 0.7 102750E— 02 0.3407863E+00 0.7102750E-02 

4 0, 7385597E+00 -0.1887922E+00 0. 7385597E*00 

5 0. 17161 12E-01 0.6304836E+00 0, 17161 12E-01 

6 -0.4901 721 E+OQ -0.2788687E+00 -0.4901721E+00 

7 -0. 1 151719E+00 -0. 6979974E-02 -0.1 1517 19E+00 

8 -0.8381367E-01 0. U26262E+01 -0.8381367E-01 

9 -0.3844098E+00 -0.2056300E-01 -0. 3844098E+00 

10 0, 1507280E+00 0.8476348E+00 0.1507280E+00 

11 ~0. 71 1 5277E+00 -0.3351524E-01 -0.7115277E+00 

12 0.3551622E+00 -0.52868B9EKK) 0. 3551622E+00 


IMAGINARY 


-0.4020393E-01 
-0.5255843E-01 
-0, 12481 90E+00 
0.1002887E-16 
-0.215612BE+00 
0. 1 227091 E^OO 
-0.5452933E-01 
0.1 5353 15E+01 
-0. 1702039E+00 
-0, 1775826E+01 
-0.2956053E+00 
0.5111248E+00 


8 

IMAGINARY 


0, 11 54731 E-01 
-0. 1781620E+00 
0. 1828979E-01 
-0. 1983724E+00 
“0. 20301 69E-16 
0. 1534676E+00 
0.1527472E+00 
0.7148006E-01 
0.5308227E+00 
-O.U17573E+00 
0. 1017405E+01 
-0.5313960E+00 


10 

IMAGINARY 


-0. 1021558E+00 
0.2439455E-17 
-0. 3407863E+00 
0. 1887922E+00 
-0.6304836E+00 
0.2788687E+00 
0.6979974E-02 
-0.1126262EKM 
0.2056300E-01 
-0.8476348E+00 
0.3351524E-01 
G.5286889E+00 
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m col 


12 


C = 1.0 (concluded) 

It 

REAL I HAG I MAR V REAL I HAG I NARY 


1 0,!85542QE-Qi 0.7251692E-Q2 0.1855420E-01 -0J251692E-02 

2 0.3952321E+00 -0. 1219082E-02 0.3952321E+0C 0.1219082E-02 

3 0.6463644E-01 0.2526831E-01 0.6463644E-01 -0.252683i£“0i 

4 0. 74259O9E+O0 ~0, 1469667E-02 0.7425909E+00 0.1469667E-02 

5 0, 1241697E+00 0.4855120E-01 0, 1241697E+00 -0.4855120E-01 

h o.moomm o.mmsE-u gjooooqoe+oi -o.3556i83E~i6 

7 0.2468634E-02 ~0,64?9524£-£>2 0.2468634E-02 0.6479524E-02 

8 -0.1603467E-02 -0.1375590E+00 -0.16O3467E-O2 0. 1375590E+00 

9 D.8601929E-02 '0.2257244E-0! 0.860I929E-02 0.2257244E~0i 

10 -0.2727015E-02 -0.2584583E+00 -0.27270I5E-02 0.25845B3E+00 

11 0,1652805E-01 -0,433627BE-0.1 0. 1652B05E-4H 0,4336278E-0I 

12 -0.2983460E-02 -0.3480553E+00 -0..29B346QE-02 0.3480553E+00 
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C * 1.50 


mm EIGENVALUES 
REAL I NAG INARY 


-0.9O992143E-O1 
-0.90992143E-O1 
-0.45491753E+OO 
-0.4549 I 753E+00 
H>. 164644 12E-01 
-0. 164644 126-0.1 
“O.B8305042E+00 
-Q.88305042E+00 
“0. 506721 02E-01 
-0/50672 102E-01 
-0.39033895E-02 
“O.39O339O3E-02 


0, 15648651E+02 
-0 . 1 564865 1E+02 
-0.61499463E+01 
0.61 4994 63E+01 
O, 17689761 E +01 
-O. 1763976 1E+01 
0. 746304 18E+00 
-0, 746304 18E+00 
-D.I1273850E+01 
0. 11273950E+01 
O. 34911 932E+00 
-0.349U932E+00 


BAHPING 

PERCENT CRITICAL 


0.58145975E+00 
0 < 581 45975E+00 
0.73769429E+01 
0. 737694 29E+01 
0. 9306 9084 E+ 00 
0. 93O69084E+0O 
0.76376673E+02 
0.76376673E+02 
0. 44901 249E+01 
0, 44901 249E+01 
0.1U79?76E+0i 
0.1U79978E+01 


SYSTEH EIGENVECTORS 


RON COL 1 

REAL IMAGINARY 


1 0. lOOOOQOE+O.l 0.8741380E-18 

2 -0. 36891 66E-05 -0.6811762E-O6 

3 -0.7146264E+00 -0.6284735E-02 

4 0.8940552E-03 0. 1759447E-03 

5 0.4323820E+00 0.9037963E-01 

6 -0.2166362E+00 -0.4527745E-01 

7 -0.9099214E-01 O.1564865E+02 

8 0.1099517E-04 -0.5766848E-04 

9 0, 1633730E+00 -0. 1118237E+02 

10 -0.2834649E-02 0.1397475E-01 

11 -0.14 53663E+01 0. 6757972E+01 

12 0. 7282431E+O0 -0.3385944E+01 


SYSTEM EIGENVECTORS 


RON COL 3 

REAL IMAGINARY 


1 -0.7OO3134E+00 -O.2558198E+0O 

2 -0.3711 115E-03 -0.1205074E-03 

3 -0.6663783E+OO -0.3852480E+00 

4 O. 13891 38E-01 0.2215325E-02 

5 O.IOOO CODE +01 -0.3258028E-16 

6 -0.5067647E+OO -0. 1049694E-02 

7 -0.1254693E+O1 0.4423267E+01 

8 -0.5722886E-O3 0.2337136E-02 

9 -0. 20661 07E+01 0.4273447E+01 

10 0.7304697E-02 -0.8643906E-01 

11 -0. 45491 75E +00 -0.6149946E+Q1 

12 0.2240806E+00 0.3M7053E+01 


REAL IMAGINARY 


O. 10OOOO0E+O1 -0.4377466E-16 
-0. 3689166E-05 0. 681 1 762E-06 
-0.7146264E+00 0.6284735E-02 
O.8940552E-03 -0.1759447E-03 
0. 4323820E+00 -0.9037963E-01 
-0.2166362E+OO 0.4527745E-01 
-0. 90992 14E-01 -0. 1564865E+02 
O. 109951 7E-G4 0.5766848E-04 
0.163373OE+OO 0, 1H8237E+02 
-0.2834649E-02 -0.1397475E-01 
-0. 1453663E+01 -O.6757972E+01 
0. 728243 1E+00 0.3385944E+01 


4 

REAL IMAGINARY 


-0. 70031 34E+00 O.2558198E+0O 
-0.3711 11 5E-03 O.12O5O74E-03 
-0. 6663783E+00 O. 3852480E+00 
0.1389138E-01 -O.2215325E-02 
O.IOO0OOOE+O1 -O.7480995E-17 
-0 . 506764 7E+00 0.1049694E-02 
-0.1254693E+01 -0.4423267E+01 
-0.5722886E-03 -0.2337136E-02 
-0, 20661 07E+01 -0.4273447E+01 
O, 730469 7E -02 O.86439O6E-01 
-0. 45491 75E+00 0.6I49946E+01 
O.224O806E+OO -0.3117053E+01 
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SYSTEM EIGENVECTORS 

C s 1.50 (continued) 

m COL 5 6 

REAL IMAGINARY REAL IMAGINARY 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


0.I703675E-01 
-0, 8833838E+00 
0.5326868E-01 
0.1Q00QQQE+01 
0.9267497E-01 
-0, 24562 16E+00 


0.3379409E-01 0.1703675E-01 -0.33794Q9E-01 
0.455776 IE-01 -O.Q83383SE+00 -0.4557761E-01 
0« 1051073E+00 0.5326868E-01 -0. 1051073E+00 
0* 26291 90E- 17 O.IGOOGQOE+Oi 0.6884684E-16 
0. 1818933E+00 0,9267497£-Ql -0.1818933E+00 
-0.1038279E+00 -0. 245621 6E+00 0. 1038279E+00 


-0.6006143E-01 O.295B121E-01 -0.6006143E-01 -0.2958121E-01 


-0, 66081 32E -01 -0.1563435E+01 -0.6608132E-01 0, 1563435E+01 
-0. i 86809 3E+00 0.9250049E-01 -0. 1868093E+00 -0.9250049E-Q1 
-0.1646441E-01 0.1768976E+01 -0.1646441E-01 -0.1768976E+01 
-0.3232907E+00 0.1609451E+00 -0/3232907E+00 -0.1609451E+00 
0, 1877131E+00 -0.4327893E+00 0.1877131E+00 0.4327893E+00 


SYSTEM EIGENVECTORS 


ROW COL 7 8 

REAL IMAGINARY REAL . IMAGINARY 


1 0. 1451437E+00 0. 1344744E-01 0, 1451437E+00 -0. 1344744E-01 
•2 -0. 3668040E-01 -0.5817737E-01 -0.366804GE-01 0.58177378-01 

3 0.5138I48E+00 G.2142O18E-01 0,513B148E+00 -0.2342018E-01 

4 -0. 158214 1E+0Q -G.8097058E-01 -0.1582141E+00 0.8097058E-01 

5 0. lOOOOOOE+Ol 0.9595189E-17 0.1000000E+01 -0.4179599E-16 

6 -0.4237223E+00 0.8672913E-01 -0.4217223E+00 -0.8672913E-01 

7 -0.1382051E+00 0.9644655E-01 -0.1382051E+00 -0.9644655E-01 

8 0.7580866E-01 0.2399881E-01 0.7580866E-01 -0,239988 LE-01 

9 -0, 46971 03E+00 0.3645470E+00 -0.4697103E+00 -0.3645470E+00 

10 Ot 2001397E+00 -0.4657473E-01 0.2001397E+00 0.4657473E-01 

11 -O.8B305O4E+OO 0.7463042E+00 -0.8830504E+00 -0.7463042E+Q0 

12 0.3076757E+00 -0.3913L93E+00 0.3076757E+00 0.3913193E+00 


SYSTEM EIGENVECTORS 


ROW COL 


10 


REAL IMAGINARY REAL IMAGINARY 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


0.5017824E-03 -0.6129145E-01 
0.1000000E+01 -0.6491661E-17 
0.2550694E-G2 -0.2044287E+00 
0.7315707E+00 O.U42539E*QO 
0.6177524E-02 -0.3781517E+00 
-0, 4778583E+00 0.1671697E+00 
-0.6912449E-01 0.2540065E-02 
-0.5067210E-01 -0.1127385E+01 
-0.2305991E+G0 0.7483220E-D2 
0.9173795E-01 -0.8305513E*00 
-0. 4266355E+00 0.1219729E-01 
0. 2126787E+00 0.5302595E+00 


0.5017824E-03 0.6129145E-01 
0.1000000E+01 0.3882799E-17 
0. 2550694E-02 0.2044287E+00 
0.7315707E+00 -0. 1142539E+00 
0.6177524E-02 0.37815 17E+00 
-0.4778583E+00 -0.1671697E+00 
-0.6912449E-01 -0, 2540065E-02 
-0. 506721 0E-01 0.U27385E+01 
-0.2305991E+00 -0.7483220E-02 
0.9173795E-01 0.8305513E+00 
-0.4266355E+00 -0.I219729E-01 
0.2126787E+00 -0.5302595E+00 
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SYSTEM EIGENVECTORS 

G = 1 .50 (concluded) 

ROW COL 11 


REAL IMAGINARY REAL 


1 0.15917G6E-01 -0.9404970E-02 0.1591786E-01 

2 0.3956620E*00 O.16O2719E-02 O.395662OE+0O 

3 O.5544906E-01 -0.3276955E-O1 0.5544906E-01 

4 0.7431O94E+00 0.1931741E-O2 O.743LO94E+0Q 

5 0, 1 065151 E+00 -0 . 62961 37E-01 O.1065151E+O0 

6 0.1O00000E+O1 -0.4163336E -16 0. 1000000E+01 

7 0.3221323E-02 0.5593942E-02 D.3221323E-02 

8 -0.21039A3E-02 O. 13ai270E+00 - 0 . 2103963E-02 

9 0.1122404E-O1 O. 1948625E-01 O.1122404E-01 

10 -0.3575O54E-02 O . 2594263E+00 -0.3575054E-02 
It O.2156526E-01 0.3743224£-0i 0.2156526E-01 
12 -O. 39O3390E-02 0.3491193E+00 -0. 3903390E-02 


12 

IMAGINARY 


0.940497OE-02 
-0.1 6027 19E-02 
0.3276955E-01 
-0.193174 IE-02 
0. 62961 37E-01 
-0.1406752E-16 
“0.5593942E-02 
-O.1381270E+0O 
-0.1948625E-O1 
-0.2594263E+0O 
-0.3743224E-01 
-0.3491193E+OO 
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C = 2.0 


SYSTEM El SEN VALUES 
REAL IMAGINARY 


-0.H801539E+00 
“0.11801 539E +00 
“0.56283407E+00 
-0. 56283407E+00 
-0.13452656E-01 
-0. 13452656E-01 
-0.17299727E+01 
"0.36949523E-01 
-0. 36949523E-0S 
-0.79872 75 1E+GQ 
-0, 439828 14E-Q2 
-0,439828i4£-O2 


0.15635587E+02 
-0. 15635587E+02 
-0,60371 4 10E+01 
0.60371410E+01 
“0. 17657839E+01 
0.17657839E+01 
0.63409056E-10 
0.11277039E+01 
-0, 1 1277039E+01 
“0.11243424E-H 
-0.35O22652E+0O 
0.35022652E+00 


DAMPING 

PERCENT CRITICAL 


0J5476556E+0O 
0.75476556E+00 
0, 92826048E+01 
0.92826048E+01 
0.76182966E+00 
0.76182966E+00 
0, 1OO0OOOOE+O3 
0,32747696E+Oi 
0, 32747696E+01 
0.10000000E+03 
0. 12557400E+01 
0. 12557400E+0I 


SYSTEM EIGENVECTORS 


RON COL 1 

REAL I HAS I NARY 


1 0. lOOOOOOE+Ol 0, 43097 04E- 17 

2 -0.3599685E-05 -0. 8S79952E-06 

3 -0.7137543E+00 “O.81679O4E-02 

4 0, 8694940E-03 0.2285858E-Q3 

5 0,4190 70 7E + 00 0, 1170330E+00 

6 “0,20996831+00 -0,5863034E“01 
1 "O.H30154E+OD 0.1563559E+02 

8 0,1 43091 4E-04 -Q.5617839E-04 

9 0.2U9440E+00 -0,ill5900Et02 

10 - 0,3 6 76 687 E - 0 2 0. 1356807E-01 

11 -0, I879336E+01 0.6538604E+01 

12 0.9414993E+00 -0.3276058EHH 


SYSTEM EIGENVECTORS 


RON COL 3 

REAL IMAGINARY 


1 -0.6837042E+00 -0.3475890E+00 

2 -0,386490 IE-03 -0. 161 6318E-03 

3 -0.6449168E+00 -0.5217492E+00 

4 0, 1428944E-01 0.2890024E-02 

5 0. lOOOOOOE+Ol 0.1764539E-16 

6 -0. 506956 2E +00 -0.1367328E-O2 

7 -0, 1713632E+01 0.4323253E+01 

8 -0.7582642E-03 0.2424267E-02 

9 -0.2786892E+01 0.41871 12E+01 

10 0, 9404902E-02 “0.8789396E-01 

11 -0,562834 1E+00 -0.6037141E+01 

12 0. 2770775E+00 0.3061336E+01 


REAL IMAGINARY 


0.1000000E+01 0.3935654E-16 
“0.3599685E-05 0.8879952E-06 
-0.7137543E+00 0.8167904E-02 
0, S694940E-03 -0.2285858E-03 
O.419O707E+OO -0. 1170330E+00 
-0. 2099683E+00 0 . 5863034E-0 I 
-0,11 801 54E+00 "0.1563559E+02 
0.1430914E-04 0,5617839E“04 
0.2U9440E+00 0.M15900E+02 
-0.3676687E-02 -0. 1356807E-01 
-0. 1879336E+01 -0.6538604E+01 
0, 9414993E+00 G.3276053E+01 


4 

REAL IMAGINARY 


-0. 6837042E+00 0. 3475890E+00 
“0, 3864901 E-03 0.1616318E-03 
-0.6449168E+00 0. 5217492E+00 
0.1428944E-01 -0.2890024E-02 
0. lOOOOOOE+Ol 0.1170938E-16 
-0. 5069562E+00 0, 1367328E-02 
-0, 1713632E+01 -0.4323253E+01 
“0, 7582642E-03 -0.2424267E-02 
-0, 27B6892E+01 -0,41871 12E+01 
0.9404902E-02 0.8789396E-01 
“0, 562834 1E+00 0,6037141E+01 
0, 2770775E+00 -0.3Q61336E+01 
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C = 2.0 (continued) 
SYSTEM EISENVECTORS 


RON COL 5 

REAL IMAGINARY REAL 


1 D. 1036344E-0.1 -0.2773020E-01 0.1036344E-01 

2 — 0. 892991 8E+00 -0.3795372E-01 -0. 892991 8E+00 

3 0.324302 IE -01 -O.8631915E-01 0.3243021E-01 

4 O.lOOOOOOE^Ol -0* 1 900064E- i 6 O.IOOOOOOE+Oi 

5 0.56467B8E-O1 -O. 1495O54E+O0 0.5646788E-01 

6 -0.2248202E+00 0.8546268E-01 -0.2248202E+00 

7 -Q.4910496E-0! -0. 1792655E-01 -0.4910496E-0i 

8 -0.5500495E-01 0.1577341E+01 -0.5500495E-01 

9 -0, 152B572E+O0 -O.56T0353E-01 -0. 1528572E+00 

10 -0.1345266E-01 -0. 1765784E+01 -O, 1345266E-01 

11 -0. 2647539E+00 -0.9769B84E-01 -0.2647539E+00 

12 0. 1539331E+00 0.3958342E+0G 0. 1539331E+00 


SYSTEH EISENVECTORS 


m COL 7 

REAL IMAGINARY REAL 


1 0.T203451E+00 -0.1651375E-T5 0.2659510E-03 

2 -0.1345592E-01 0.26978OOE-17 O.IOOOOOOE+Ol 

3 0.4739241E+00 KM695828E-15 O. 1352447E-02 

4 -O.9214683E-01 O, 1402348E-16 0.7296491E+00 

5 0.1000000E+01 -0. 1637145E-16 0.3276075E-02 

6 -O.4416152E+O0 -0.3260738E-16 -Q.4745571E+00 

7 -0.2081937E+Q0 -0.4868047E-14 -O.5G3O5O6E-01 

8 0.3192824E-O1 -0.5676815E-17 -0.3694952E-01 

9 “O.8I98757E+0O 0.2834836E-14 -0. 167794 1 E+00 

10 0. 1594115E+0O -0. 1818072E-16 0. 6701S43E-01 

11 -0. 1729973E+01 -0. 1562552E-14 -0,3104002E+0G 

12 0. 7639823E+O0 0.8151574E-15 0. 1546775E+00 


SYSTEM EISENVECTORS 


RON COL 9 

REAL IMAGINARY REAL 


1 0. 2659510E-03 -0.4459968E-01 0, 1416727E+00 

2 O.10O0OO0E+01 -0.96765O4E-17 -0, 5840039E-01 

3 D , 1 352447E-02 -0 . 1487484E+00 0. 5081991E+00 

4 D.7296491E+O0 0.8333625E-O1 -0.154Q582E+00 

5 0, 3276075E-02 -0, 2751424E+00 0. IDOOOOOE+Ol 

6 -0.47 4557 1E+00 0. 1216124E+00 -O.3479999E+0O 

7 -0.5030506E-01 0.1348023E-O2 -O, 1131579E+0D 

8 -0. 3694952E-01 -8.1127704E+O1 0, 4664600E-01 

9 -0. 167794 1E+0O 0.3971022E-02 -G. 40591 26E+O0 

10 0.6701S43E-01 -O.8259074E+O0 0. 1230505E+00 

11 -0.31O4O02E+OO 0.6471936E-02 -0. 79S7275E+00 

12 0. 1546775E+00 0.53066641+00 0, 2779571E+00 


6 

IMAGINARY 


0.2773020E-O1 
0.3795372E-01 
0.8631 9 15E-01 
-0.4100995E-16 
O. 1495O54E+O0 
-0.8546268E-01 
0, 1792655E-01 
-0.1577341E+01 
0.5610353E-01 
0.1765784E+01 
0.97698S4E-O1 
-0.3958342E+O0 


8 

IMAGINARY 


0.4459968E-01 
0.9690057E-i8 
0.1487484E+OO 
-0.8333625E-01 
0.2751424E+0G 
-O. 12161 24E+O0 
-O.1348023E-02 
0. 1127704E+01 
HU971022E-02 
O. 8259O74E+0O 
-0.6471936E-O2 
- 0 . 5 3 0 6 6 6 4 E +00 


10 

IMAGINARY ' 


0.73O4812E-17 
0.4336809E-17 
-0.45S3465E-I6 
0 . 474 33 8 5£ - 1 8 
0.1897354E-18 
-0.5468445E-16 
-0.3152938E-14 
0.238931 IE- 16 
0.1725196E-14 
O.2708473E-16 
-0.9 10296 IE- 15 
0.4483447E-15 
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C = 2.0 (concluded) 
SYSTEM EIGENVECTORS 


ROW COL 11 

REAL IMAGINARY REAL 


1 0, 1322368E-01 0.1050203E-01 0.I322368E-03 

2 O.3961138E+O0 -0. 18I5051E-02 0,3961 138E+00 

3 0,46061 I9E-01 0.3658992E-01 0,4606119E-01 

4 0, 7436541£-*'00 -0.2187163E-02 0. 7436541E+00 

5 0.8847686E-01 0.7029820E-01 0.8847686E-01 

6 0, 1000000E+01 0.2986977E-U 0.1000000E+01 

7 0.36W30E-02 -0. 4677474E-02 0. 3619930E-02 

8 -0.2377899E-02 -0. 1387216E+00 -0,2377899E-02 

9 0,12612 17E -0 1 -0 , 1 6292 78E-0 1 0, 126121 7E-01 
tO -0,4O36803E-02 -0.2604378E+00 -0.4036803E-02 

11 0.24231 15E-01 -O.3129613E-01 0. 24231 15E-01 

12 -O. 4398281 E *02 -0.3502265E+OO -0.4398281E-O2 


12 

IMAGINARY 


-O. 1O5O203E-01 
0.1815051E-02 
-0.3658992E-O1 
0.2187163E-02 
-0. 7029820E-01 
0.5B27587E-I8 
0.4677474E-02 
0,1 3872 16E+O0 
0. 1629278E-01 
0, 2604378E+00 
O.3129613E-01 
0, 3502265E+00 
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